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APPENDIX  D 


Emissions  Data 
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E  -  1 
E.l   THE  REQUIREMENT 


The  emissions' ^<3ata  supplied  by  Senes  Consultants  Ltd.  (see 
Appendix  D)  include  two  broad  categories  of  contaminants,  namely 
particulate  matter  (PM)  and  volatile  organic  compounds  (VOC's). 
These  broad  categories  include  dozens  (PM)  or  hundreds  (VOC)  of 
different  chemicals. 

Reduced  emissions  of  particulate  matter,  regardless  of  the 
composition,  can  yield  health  (respiratory  problems,  eye 
irritations  etc.)  and  other  (soiling,  visibility  reduction) 
benefits.   But  most  health  and  environmental  benefits  derived 
from  reduced  emissions  of  PM's  and  VOC's  result  from  lower 
concentrations  of  specific  chemicals.   Thus,  the  composition  of 
the  VOC  and  PM  emissions  must  be  determined  so  that  the  changes 
in  emissions  of  specific  chemicals  can  be  examined.   This  was  one 
of  the  major  challenges  of  the  project. 


E . 2   APPROACH 


Not  surprisingly,  the  composition  of  the  VOC's  (PM's)  emitted 
varies  with  the  specific  industrial  process.   The  emissions  from 
a  coal-fired  boiler,  a  gas-fired  boiler  and  an  oil-fired  boiler 
are  different.   To  determine  the  composition  of  VOC  (PM) 
emissions,  it  is  necessary  to  begin  with  the  processes  used  in 
each  economic  sector. 

E.2.1   Speciation 

The  only  source  of  detailed  data  on  the  composition  of  VOC's  and 
PM's  is  of  the  United  States  Environmental  Protection  Agency's 
(EPA) ,  National  Emissions  Data  System  (NEDS) . 
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As  in  Canada,  U.S.  industries  are  classified  by  Standard 
Industrial  Classification  (SIC)  codes.   Unfortunately,  U.S.  SIC 
codes  differ' f torn  Canadian  SIC  codes. 

In  the  NEDS  database  each  U.S.  SIC  is  divided  into  processes 
which  are  regarded  as  sources  of  emissions  and  are  assigned  a 
unique  Source  Classification  Code  (SCC) .   To  use  the  SCC 
information  on  PM  and  VOC  composition  we  needed  to: 

.  select  the  SCC's  that  correspond  to  each  of  the  Canadian 
industries  or  industry  segments  covered  by  this  study; 

.  ensure  that  the  SCC's  selected  correspond  to  the  processes 
and  emissions  estimated  by  Senes;   and 

.  aggregate  the  SCC  profiles  of  PM's  and  VOC's  into  profiles 
corresponding  to  the  Canadian  SIC  codes  covered  by  this 
study . 

These  steps  are  described  below. 


E.3   IDENTIFICATION  OF  SCC'S 


The  first  task  is  to  identify  the  unique  Source  Classification 
Codes  (SCC's)  needed  for  our  analysis.   The  Canadian  SIC  codes 
for  the  industries  or  industry  segments  covered  by  the  study  are 
known.   Statistics  Canada  provides  a  conversion  table  that 
relates  each  Canadian  SIC  code  to  the  corresponding  U.S.  SIC 
codes.   We  used  this  conversion  table  to  identify  the  relevant 
U.S.  SIC  codes.   The  number  of  U.S.  SIC  codes  identified  was 
approximately  50  per  cent  larger  than  the  corresponding  number  of 
Canadian  SIC  codes. 
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We  then  obtained  from  the  National  Emissions  Data  System  (NEDS) 
data  for  the  SCO's  corresponding  to  the  identified  U.S.  SIC 
codes.   We  rec'eived.  data  on  approximately  15,000  SCC's. 


E.4   sec  DATA 


The  data  obtained  for  each  SCC  included: 

.  the  U.S.  SIC  code  that  uses  the  process; 

.  the  type  of  process; 

.  a  description  of  the  process  including  notes  such  as  the 

fuel  used; 
.  the  corresponding  speciation  profile  number  which  identifies 

the  specific  chemicals  that  comprise  the  VOC  (PM)  emissions 

for  that  process;   and 
.  the  total  VOC  (PM)  emissions  from  the  given  process  by  all 

plants  in  the  United  States  in  1985. 

In  theory  each  SCC  has  a  unique  chemical  breakdown  for  its  VOC 
emissions  and  its  PM  emissions.   In  practice,  generalized 
profiles  are  common  to  many  SCC's.   For  example,  the  composition 
of  the  VOC  and  PM  emissions  from  oil-fired  boilers  might  be 
similar  for  many  industries  (SIC's)  that  use  such  boilers. 

It  is  also  important  to  understand  that  the  SCC's  corresponding 
to  a  given  SIC  code  do  not  represent  a  single  set  of  processes 
found  in  all  plants  in  that  industry.   Sticking  with  the  boiler 
example,  a  given  SIC  code  may  have  SCC's  that  correspond  to 
oil-fired,  coal-fired,  and  gas-fired  boilers.   A  given  plant 
might  use  only  a  single  fuel  for  its  boiler,  so  only  one  of  the 
SCC's  would  be  relevant  for  its  emissions. 
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E.5   ELIMINATION  OF  MINOR  SCC'S 


As  mentioned  above,  a  large  number  (more  than  15,000)  processes 
(SCC's)  are  sources  of  emissions  in  the  industrial  sectors  we 
considered.   We  had  to  find  some  way  of  cutting  down  this  number 
to  make  the  calibration  of  SCC's  with  the  industry  sectors 
covered  by  the  study  and  the  emissions  data  provided  by  Senes  a 
manageable  task. 

We  used  the  following  procedure.   We  eliminated  all  SCC's  whose 
contribution  to  total  PM  or  total  VOC  emissions  in  the  (U.S.) 
SIC  in  question  was  less  than  5  percent.^   In  most  cases,  we  were 
left  with  a  small  group  of  SCC's  which  contributed  75  percent  or 
more  of  total  PM  or  VOC  emissions  of  U.S.  firms  in  that  industry. 

In  a  few  instances  we  were  left  with  less  than  75  percent  of 
total  PM  or  VOC  emissions  using  the  "less  than  5  percent"  rule. 
In  those  cases  we  conducted  an  ordered  search  for  "minor"  SCC's 
and  put  them  back  in  the  SCC  grouping  for  the  American  SIC  in 
question.   We  kept  adding  "minor"  SCC's  until  we  had  75  percent 
of  total  PM  and  total  VOC  emissions  accounted  for.   Our  ordered 
search  involved  looking  for  the  most  important  SCC's  of  less  than 
5  percent,  i.e.,  4-5  percent,  the  3-4  percent,  etc. 

The  process  of  eliminating  minor  SCC's  reduced  the  number  from 
over  15,000  to  less  than  1,000. 


^   Only  if  the  SCC  contributed  less  than  5  percent  of  both 
total  VOC  and  PM  emissions  was  it  excluded.   If  the  SCC 
contributed  5.1  percent  of  total  VOC  emissions  and 
1.0  percent  of  total  PM  emissions  it  was  retained. 
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E.6   ELIMINATION  OF  SUB-SECTORS  NOT  COVERED 
BY  THIS  STUDY 


The  American  SIC  codes  include  sub-sectors  that  are  not  found  in 
Ontario  or  were  excluded  from  this  study.   For  example,  the 
American  SIC  code  for  fish  processing  includes  canneries  which 
are  not  found  in  Ontario.   Hence,  the  canneries  sub-sector  and 
related  SCC's  were  excluded.   Similarly,  only  the  coating 
operations  of  the  automotive  and  automotive  parts  industries  are 
addressed  in  this  study.   Hence,  SCC's  relating  to  other 
processes  in  the  automotive  and  automotive  parts  industries  were 
eliminated. 

Few  problems  were  encountered  in  matching  the  sector  and 
sub-sector  definitions  to  those  within  the  scope  of  this  study. 
Matching  the  process  definitions,  with  the  exception  of  boilers, 
also  did  not  create  any  problems.   Where  questions  arose  we 
contacted  Senes  Consultants  to  ensure  that  the  sub-sector  and/or 
process  coverage  in  terms  of  SCC's  matched  emission  sources  they 
identified  in  Ontario. 


E.7   BOILERS 


Boilers  are  an  integral  part  of  many  industrial  processes.   They 
are  also  used  to  provide  space  heating  for  larger  buildings.   It 
is  not  surprising  then  to  find  that  a  number  of  the  SCC's  related 
to  boilers.   However,  Senes  advised  us  that  boilers  are  not 
important  sources  of  emissions  in  many  of  the  sectors  studied  and 
hence  boiler  emissions  were  not  included  in  their  estimated 
emissions  for  these  sectors. 
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Boilers  therefore  presented  three  problems  for  the  analysis: 

.  Whether,  for  a  particular  economic  sector,  the  boiler 
aspect,  Should  be  included  or  not. 

.  Whether  boilers  used  for  space  heat  should  be  considered. 

.  If  it  was  decided  to  include  boiler  emissions  in  a 
particular  sector,  how  the  speciation  of  VOC  and  PM 
emissions  obtained  from  the  U.S.  should  be  adjusted  to 
account  for  the  greater  use  of  natural  gas  and  oil  in 
Ontario  relative  to  the  United  States. 

Our  approach  to  these  problems  is  described  below.   The  approach 
is  apparently  simplistic,  but  we  considered  it  better  to  be 
consistent  in  a  project  with  uncertainty,  in  any  case,  so 
prevalent. 

E.7.1   Boiler  Inclusion  in  the  Analysis 

The  EPA  NEDS  database  may  indicate  that  boilers  are  sources  of 
VOC's  and  PM's  in  a  particular  U.S.  industrial  sector. 
Regardless  of  that  information,  we  followed  the  advice  of  the 
consultants  who  provided  the  emissions  data  (Senes)  as  to  whether 
boilers  had  in  fact  been  included  in  their  calculations.   If 
boiler  emissions  were  included  in  the  emissions  estimates,  EPA 
boiler  speciations  were  used  for  PM's  and  VOC's  emitted  by  the 
sector  in  question.   If  boiler  emissions  were  not  included  in  the 
emissions  estimates,  the  SCO's  relating  to  boilers  were  dropped. 

Eliminating  the  emissions  from  boilers  had  the  effect  of 
reducing  the  percentage  of  total  VOC  or  PM  emissions  covered  by 
the  selected  SCC's  below  75  percent  of  total  VOC  and  PM  emissions 
by  that  industry  in  most  cases. 
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E.7.2   Boiler  Use  for  Space  Heat 

Emissions  from  boilers  used  for  space  heating  were  excluded  in 
all  cases ,  ekbept  to  the  extent  that  the  boiler  would  be  included 
as  part  of  the  general  boiler  use  in  the  sector. 

E.7.3   Boiler  Fuel  Conversion,  Sectors 

We  judged  it  necessary  to  make  boiler  fuel  conversions  in  three 
sectors:   electric  power  generation,  pulp  and  paper,  and 
petroleum  refineries. 

Ontario  Hydro  information  for  the  six  power  generation  plants 
examined  in  this  study  indicates  that  only  three  fuels  are  used: 
lignite  coal,  bituminous  coal,  and  heavy  fuel  oil.   The  EPA  data 
also  indicated  boilers  fueled  by  wood  bark  waste,  distillate  oil 
and  natural  gas.   We  treated  these  fuels  as  heavy  fuel  oil  and 
adjusted  emissions  as  described  below. 

In  the  pulp  and  paper  sector,  we  assumed  that  wood  bark  waste, 
natural  gas,  and  heavy  fuel  oil  are  used  in  Ontario.   Where  EPA 
identified  the  use  of  coal,  we  treated  these  emissions  as  if 
they  came  from  heavy  fuel  oil.   For  petroleum  refineries,  we 
assumed  natural  gas  and  heavy  fuel  oil  are  used  to  fuel  boilers. 
Where  EPA  referred  to  use  of  light  oil,  we  converted  this  to 
heavy  fuel  oil. 

E.7.4   Boiler  Fuel  Conversion,  Procedure 

Next  we  had  to  convert  emissions  from  one  fuel  type  to  another. 
We  used  the  following  procedure. 

Senfes  provided  us  with  information  on  VOC  and  PM  emissions  per 
tonne  of  steam  for  different  fuels  by  boiler  type  (i.e.,  the  many 
different  kinds  of  boiler  that  are  powered  by  one  type  of  fuel) . 
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VOC  and  PM  emissions  per  tonne  of  steam  do  not  vary  greatly  by 
boiler  type.   The  emissions  per  tonne  of  steam  depend  much  more 
on  the  fuel  used  than  on  the  boiler  design.   We  calculated 
typical  VOC  and  PM  emissions  per  tonne  of  steam  by  fuel  type, 
regardless  of  the  boiler  type. 

To  convert  the  EPA  data  from  one  boiler  fuel  to  another,  we 
started  from  the  premise  that  the  steam  required  would  not 
change.   Then  the  quantity  and  composition  of  the  VOC  and  PM 
emissions  resulting  from  the  generation  of  that  steam  were 
adjusted  to  reflect  the  different  boiler  fuel. 

For  example,  assume  the  NEDS  data  for  the  pulp  and  paper  industry 
indicated  that  particulate  emissions  from  bituminous  coal-fired 
boilers  amounted  to  500  tonnes  per  year.   The  steam  generated  was 
estimated  at  (500  x  1,000/5.35  =)  93,458  tonnes  per  year.2 
Assuming  this  same  93,458  tonnes  of  steam  was  generated  by  a 
heavy  oil  fired  boiler,  the  PM  emissions  would  be  (93,458 
X  0.075/1,000  =)  7  tonnes  per  year.-^  The  magnitude  of  this  shift 
in  PM  emissions  —  from  500  tonnes  to  7  tonnes  in  the  example  — 
indicates  why  we  felt  it  important  to  adjust  the  mix  of  boiler 
fuels  where  possible. 

The  PM  emissions  estimated  by  converting  from  one  boiler  fuel  to 
another  (bituminous  coal  to  heavy  oil)  were  added  to  the 
emissions  reported  by  the  NEDS  database  for  the  substitute  fuel 
(heavy  oil) .   The  SCC  profile  for  boilers  using  that  fuel  was 
assumed  to  apply  to  the  revised  total  emissions. 


2  500  tonnes  of  PM  x  1,000  kg/tonne  divided  by  PM  emissions 
of  5.35  kg  per  tonne  of  steam. 

^    93,458  tonnes  of  steam  x  0.075  kg  of  PM  emissions  per 
tonne  of  steam/1,000  kg  per  tonne  =  7  tonnes  of  PM 
emissions. 
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These  calculations  eliminate  some  SCC  profiles  —  those  for  the 
boiler  fuels  that  have  been  converted  —  from  the  set  of  profiles 
for  the  industry  and  assign  higher  total  PM  and  VOC  emissions  to 
other  SCC  profiles  : —  those  for  the  substitute  fuels  —  for  that 
industry.   The  total  emissions  of  PM  and  VOC  by  SCC  obtained  from 
the  NEDS  database,  adjusted  for  the  different  mix  of  boiler 
fuels,  are  used  as  weights  for  the  SCC  profiles  in  the 
calculation  of  overall  PM  and  VOC  emissions  profiles  for  the 
industry. 

This  process  may  appear  unsatisfactory.   However,  the  reader 
should  remember  that  the  objective  of  the  speciation  process  is 
to  define  the  characteristic  proportions  of  contaminants  in  the 
emissions  from  a  particular  economic  sector.   Thus,  apparent 
errors  resulting  from  our  adjustments  may  alter  the  composition 
of  the  VOC  or  PM  emissions  speciations  somewhat,  but  will  not 
affect  the  total  amount  of  such  emissions.   The  amounts  are  those 
estimated  by  Senes. 


E.8   CALCULATION  OF  SPECIATIONS 


Having  identified  the  SCO's  that  account  for  the  bulk  of  the  VOC 
and  PM  emissions  in  the  sectors  and  sub-sectors  of  interest  for 
this  study  and  having  adjusted  the  process  fuels  where 
appropriate,  it  remained  to  calculate  a  VOC  and  PM  speciation 
for  each  industry.   That  was  done  as  follows. 

All  of  the  selected  SCO's  corresponding  to  a  given  Canadian  SIC 
were  identified.   The  composition  of  each  of  these  SCO's  was 
determined  from  the  NEDS  database.   This  lists  every  chemical 
that  is  found  in  the  VOC  (PM)  emissions  by  that  SCC  together  with 
the  percent  weight  of  that  chemical  of  total  VOC  (PM)  emitted  by 
that  process  (SCC)  by  U.S.  plants  in  1985,   Next,  weighted 
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averages  of  the  same  contaminant  from  different  SCO's  were 
combined.   This  total  was  used  to  determine  each  contaminant's 
percentage  sha're  of  the  total  VOC  (PM)  emissions  for  a  specific 
Canadian  SIC'. 

In  a  few  instances,  the  emissions  data  provided  by  Senes 
Consultants  (Appendix  D)  identifies  separately  chemicals  that 
are  constituents  of  the  PM  or  VOC  emissions.   We  checked  with 
them  to  determine  whether  those  chemicals  were  emitted 
separately,  in  addition  to  the  emissions  that  of  that  chemical  as 
part  of  the  PM  or  VOC  emissions.   Where  the  emissions  were  said 
to  be  part  of  the  PM  or  VOC  emissions,  the  quantities  of  that 
chemical  estimated  by  the  speciation  of  the  PM  or  VOC  emissions 
were  subsequently  replaced  by  Senes'  estimate  of  emissions  for 
the  contaminant.   This  uses  the  more  relevant  information 
generated  by  Senes,  and  eliminates  double  counting  of  emissions 
of  those  contaminants. 


E .  9   RESULTS 


The  PM  and  VOC  speciations  for  the  Canadian  SIC  codes  used  in 
this  analysis  are  shown  below. 

The  profiles  break  down  the  PM  and  VOC  emissions  into  their 
constituent  chemicals.   Some  of  those  constituent  chemicals  are 
Regulation  3  08  Contaminants.   Others  are  not  on  the  list  of 
96  chemicals  we  were  asked  to  consider.   Each  profile  shows  the 
percentage  contribution  of  each  Regulation  3  08  Contaminant  and  a 
balance  percentage  of  "other"  chemicals  making  up  the  profile. 

This  "other"  category  typically  accounts  for  most  —  usually  over 
60  percent  —  of  the  PM  or  VOC  emissions  in  a  sector.   The 
relatively  large  size  of  this  "other"  category  raised  the  concern 
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that  potentially  harmful  contaminants  were  being  overlooked.   To 
allay  this  concern  we  provided  a  list  of  all  VOC  and  PM  chemicals 
on  the  NEDS  dcitabase  to  the  Ministry  of  the  Environment  for 
review.   The' -Ministry  was  satisfied  that  the  list  of  contaminants 
originally  specified  identified  the  chemicals  of  concern.   This 
review  led  to  the  addition  of  three  chemicals  to  the  list  of 
contaminants . 

Thus,  the  "other"  category  can  be  considered  to  consist  of 
relatively  harmless  chemicals.   The  steps  taken  to  reduce  the 
more  hazardous  chemicals  in  the  VOC  or  PM  emission  streams  will 
have  the  effect  of  lowering  emissions  of  these  more  benign 
constituents  as  well  at  no  additional  cost.   Because  of  the 
relatively  benign  nature  of  these  "other"  chemicals  no  benefits 
are  deemed  to  be  derived  from  their  curtailment. 
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OEEN  PIT  MINING 
SIC  0611  Gold  Mines 

m.  Contaminant  ;,iv  Per  Cent 


SILICA 

OTHER 

Subtotal 


31.42 

68.58 

100.00 


SIC  0617  Iron  Ore  Mines 


EM  Contaminant 

Per  Cent 

ARSENIC 

2.64 

CAEMIUM 

0.75 

CHRCMIUM 

0.23 

FLUORIDE 

0.52 

MERCURY 

0.01 

MANGANESE 

0.68 

NICKEL 

0.14 

LEAD 

4.93 

SFTJTCTUM 

0.05 

Sn.TCA 

4.54 

TIN 

0.20 

OTHER 

85.30 

Subtotal 

99.99 

POCO  INDDSTRTES 

SIC  lOU  Meat  and  Meat  Products 

IM  Contaminant    Per  Cent 


SILICA 

OTHER 

Subtotal 


15.00 

85.00 

100.00 


SIC  1012  Poultry  Products 

IM  Contaminant    Per  Cent 


SILICA 

OTHER 

Subtotal 


15.00 

85.00 

100.00 


VOC  Contaminant 

BENZENE 

ETHYLBENZENE 

ETHYLENE 

PORMAIDEHYDE 

TOLUENE 

OTHER 

Subtotal 


VDC  Ccrrtaminant 

OTHER 
Subtotal 


Per  Cent 

0.11 
0.01 
0.63 
0.81 
0.04 
98.40 
100.00 


Per  Cent 

100.00 
100.00 
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SIC  1021    Fish  PixxJucts 

m.  Contaminant    Per  Cent 


SniCA 

OTHER 

Subtotal 


15.00 
■  85.00 
100.00 


SIC  1049  Other  Dairy  Products  (Whey 
Ecying) 


FM  Contaminant 

Per  Cent 

ARSENIC 

0.23 

CAEMIUM 

0.07 

CHPCMIUM 

0.03 

FUUORIDE 

0.05 

MERCURY 

0.01 

MANGANESE 

0.07 

NICKEL 

0.03 

LEAD 

0.68 

SFJ  FKHM 

0.01 

SILICA 

13.76 

TIN 

0.03 

OTHER 

85.04 

Subtotal  1 

100.01 

RJLP  AND  PAHilR  AND  AI.I.IKI)  P 

VDC  Contaminant 

OTHER 
Subtotal 


SIC  2711  Pulp  Industry 

IM  Contaminant    Per  Cent 


ARSENIC 

0.01 

CAEMIUM 

0.01 

CHRCMIUM 

0.03 

FUJORIDE 

1.19 

MEROJRY 

0.08 

MANGANESE 

0.15 

NICKEL 

0.03 

LEAD 

0.01 

SKTKNIUM 

0.02 

SILICA 

0.55 

TIN 

0.02 

OTHER 

97.90 

SiJbtotal 

100.00 

VOC  Contaminant 

BENZENE 

EIHYLBENZENE 

ETHYLENE 

FORMALDEHYDE 

TOLUENE 

OTHER 

Subtotal 


Per 

Cent 

100. 

,00 

100, 

,00 

Per  Cent 


62 
48 
67 
48 


1.00 
74.74 
99.99 
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SIC  2712     Newsparirrt  Industry 
EM  Contaitiinant  Per  Cent 


ARSINIC 

'  ';.i,  0.01 

CACMIUM 

-■■■  -  ■  0.01 

CHRCMIUM 

0.02 

FIDORIDE 

0.85 

MERCURY 

0.05 

MANGANESE 

0.12 

NICKEL 

0.43 

LEAD 

0.02 

SETFNIUM 

0.01 

STT.TCA 

2.03 

TIN 

0.02 

OTHER 

96.43 

Subtoti^l 

100.00 

VOC  Contaminant 

Per  Cent 

1,1/  l^TRICHIDRDEIHANE 

0.02 

1,3-BOTADIENE 

0.06 

acryldnitkhe 

0.05 

BENZENE 

9.46 

BUTYL  ACRYIATE 

0.03 

BUTYL  CKTJOSOLVE 

0.02 

CARBON  I'h'i'RACHDDRIDE 

0.04 

CKTiaSOLVE 

0.02 

rrr.TncrjTVF  ArrrrATF: 

0.02 

CHLOROFORM 

0.03 

EPICHDDRDHYDRIN 

0.02 

ETHYT BENZENE 

0.73 

ETHYLENE 

7.48 

ETHYLENE  DIBRCMIDE 

0.02 

ETHYLENE  DICHIDRIDE 

0.09 

ETHYLENE  OXIDE 

0.02 

FORMALDEHYDE 

4.55 

Mh'i'HYL  ETHYL  KETONE 

0.11 

Mh'i'HYL  ISOBUTYL  KETONE 

0.04 

NAPTTTATENE 

0.03 

PKHCHJJJi<OETEiYIENE 

0.04 

PHENOL 

0.55 

EROPIONAIDEHYDE 

0.02 

PROPIONIC  ACID 

0.02 

PROPYLENE  OXIDE 

0.03 

STYRENE 

0.17 

TOLUENE 

2.16 

TOLUENE  DIISOCYANATE 

0.04 

TRICHLDROETHYLENE 

0.02 

VINYL  CHLORIDE 

0.06 

OTHER 

74.06 

Subtotu^l 

100.01 
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SIC  2713  P^jerixard  Inctustiy 

RI  Contaminant    Per  Cent 


ARSENIC 

;.l>'  .  0.01 

CAEMIUM 

■  -  '  o.oi 

CHROMIUM 

0.01 

FIUORIDE 

1.18 

MERCURY 

0.08 

MANGANESE 

0.17 

NICKEL 

0.50 

LEAD 

0.02 

SILICA 

2.06 

TIN 

0.02 

OTHER 

95.92 

Subtotal  1 

99.98 

SIC  2811  Business  Forms  Printing 


SIC  2819  Other  Comnercial  Printing 

FM  Contaminant    Per  Cent 


ARSENIC 

1.04 

CAEMIUM 

0.34 

CHROMIUM 

0.12 

FUJORIDE 

0.22 

MERCURY 

0.05 

MANGANESE 

0.32 

NICKEL 

0.12 

LEAD 

3.10 

SETiENIUM 

0.03 

SILICA 

8.36 

TIN 

0.15 

OTHER 

86.17 

Subtotal 

100.02 

VOC  Contaminant 

BENZENE 

ETHYLBENZENE 

ETHYLENE 

PORMAIDEHYDE 

TOLUENE 

OTHER 

Svibtotal 


Per  Cent 

11.08 

0.43 

9.75 

6.69 

0.76 
71.27 
99.98 


VOC  Contaminant 

Per  Cent 

ML'i'HYL  EIHYL  KETONE 

54.27 

MEl'HYL  ISOBUTYL  KETONE 

14.78 

TOLUENE 

6.48 

OTHER 

24.47 

Subtotal 

100.00 

VOC  Contaminant 

Per  Cent 

BENZENE 

0.04 

BUTYL  rmnsoLVE 

0.11 

CKTJOSOLVE 

0.01 

ETHYT, BENZENE 

3.77 

ETHYLENE 

0.12 

ETHYLENE  DIBRCMIDE 

0.07 

FORMALDEHYDE 

1.22 

MLTHYL  ETHYL  KETONE 

5.51 

ML'i'HYL  ISOBUTYL  KETONE 

2.71 

TOLUENE 

3.73 

OTHER 

82.67 

Subtotal 

99.96 
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SIC  2821  Plate  Making,  Typesetting  and 
Bindery 


FM  Contaminant 

J  Per  Cent 

CAEMIUM 

V   o.Ql 

FLUORIDE 

0.36 

MANGANESE 

0.04 

NICKEL 

0.01 

lEAD 

0.01 

SHIGA 

0.21 

TIN 

0.01 

OTHER 

99.35 

S\abtc±al 

100.00 

VDC  Contaminant 

Per  Cent 

1,1, 1-TRICHIORDETHANE 

0.03 

1,3-BUIADIENE 

0.05 

ACRyioNi'rkTTF; 

0.04 

BENZKNE 

0.20 

BUTYL  ACRYIAIE 

0.03 

BUTYL  CKTIOSOLVE 

6.06 

CARBON  'i'LTKACHLDRIDE 

0.05 

CETinSOTVE 

0.03 

rPT,Tr«nTVF  acfttatr 

0.03 

CHLOROFORM 

0.04 

EPICHLORDHYDRIN 

0.03 

EIHYT,RFNZi:2re: 

0.55 

ETHYLENE 

0.20 

ETHYLENE  DIHRCMIDE 

0.03 

ETHYLENE  DICHLORIDE 

0.04 

ETHYLENE  OXIDE 

0.03 

FORMALDEHYDE 

0.11 

Mh'lHYL  ETHYL  KETONE 

0.59 

MhTHlsIL  ISUBLTTYL  KETONE 

0.37 

NAFTTTATITNE 

0.01 

PERCHLORDETHYLENE 

0.05 

HIENOL 

0.03 

PROPIONALDEHYDE 

0.03 

PROPIONIC  ACID 

0.03 

PROPYLENE  OXIDE 

0.04 

STYRENE 

0.08 

TOLUENE 

35.36 

TOLUENE  DIISOCYANATE 

0.05 

TRICHLOROETHYLENE 

0.03 

VINYL  CHLORIDE 

0.03 

OTHER 

55.77 

Subtotal 

100.02 

SIC  2841  Newspaper,  ffegazine  and 
Periodical  Publishing 


VDC  Contaminant 

Per  Cent 

BUTYL  CFTTOSOLVE 

0.43 

EIHYTRKNZENE 

0.28 

ETHYLENE 

0.03 

ETHYLENE  DIBRCMIDE 

0.14 

FORMALDEHYDE 

0.36 

ML'IHYL  ETHYL  KETONE 

12.16 

ML'iHYL  ISOBUTYL  KETONE 

4.79 

TOLUENE 

6.50 

OTHER 

75.31 

Subtotal 

100.00 
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SIC  2849  Other  Ocnfcined  Publishing  and 
Printing  Industries 


VOC  Contaminant 
EIHYLBENZENE 
OTHER 
Subtotal 


Per  Cent 

6.00 

94.00 

100.00 


IRCN  AND  •'^•'I'HIKI.  INDUSTRIES 
SIC  2911  Ferro-Alloys 

EM  Contaminant    Per  Cent 


ARSENIC 

0.21 

CAEMTUM 

0.20 

CHRCMIUM 

1.38 

MANGANESE 

0.29 

NICKEL 

0.10 

LEAD 

3.56 

SELENIUM 

0.02 

SILICA 

6.61 

TIN 

0.01 

OTHER 

87.61 

Subtotal 

99.99 

VOC  Contaminant 

Per  Cer 

1,1, 1-IRICHDDRDETHANE 

0.04 

1,3-BUIADIENE 

0.38 

ACRYLDNl'iRIIE 

0.05 

BENZENE 

9.28 

BUTYL  ACRYIATE 

0.04 

BUTYL  nTTTOSOLVE 

0.04 

CARBON  'i'Ei'RACHLDRIDE 

0.06 

rvTJnPDJ^TF. 

0.03 

(r.\  ,1  («OLVE  ACh'iATE 

0.03 

CHIORDFORM 

0.05 

EPICHIDRDHYDRIN 

0.04 

ETHYTiRFNZENE 

0.06 

ETHYLENE 

19.56 

ETHYLENE  DIBRCMIDE 

0.04 

ETHYLENE  DICHLDRIDE 

0.05 

ETHYLENE  OXIDE 

0.04 

FORMAIDEHYDE 

0.15 

ML'iHYL  ETHYL  KETONE 

0.11 

Mh'i'HYL  ISOBUTYL  KETONE 

0.05 

NAPTHATFMP 

0.02 

PEKCHUJROETHYLENE 

0.07 

PHENOL 

0.04 

PRDPIONAIDEHYDE 

0.04 

PROPIONIC  ACID 

0.04 

PROPYLENE  OXIDE 

0.05 

STYRENE 

0.11 

TOUJENE 

0.64 

TOLUENE  DIISOCYANATE 

0.06 

TRICHIDROETHYLENE 

0.04 

VINYL  CHLORIDE 

0.04 

OTHER 

68.75 

Subtotal 

100.00 

E  -  18 


SIC  2912  Steel  Foundries 


HI  Contaminant 

Per  Cei 

ARSENIC     •  ';., 

l",   0.02 

C?^EMIUM     -••  - 

0.03 

CHRCMTUM 

1.36 

MANGANESE 

3.35 

NICKEL 

1.04 

LEAD 

0.50 

snicA 

2.92 

TIN 

0.01 

OTHER 

90.75 

Subtotal 

99.98 

VDC  Contaminant 

Per  Cent 

1,1,  l-TRICHDDROETHANE 

0.43 

1,3-BLnADIENE 

0.67 

ACRYIDNi'IklLE 

0.58 

BENZENE 

2.82 

BUTYL  ACRYIATE 

0.43 

BUTYL  CELDDSOLVE 

0.44 

CARBON  TLTRACHTORIDE 

0.65 

CFTJaSOTVK 

0.36 

fKI.KiWm.VH  ACk'I'ATR 

0.36 

OffORDFORM 

0.51 

EPICHLDRDHYDRIN 

0.42 

EIHYT BENZENE 

0.65 

ETHYLENE 

2.86 

ETHYLENE  DIHHCMIDE 

0.42 

ETHYLENE  DICHIORIDE 

0.58 

ETHYLENE  OXIDE 

0.38 

FORMALDEHYDE 

1.55 

ML'IHYL  ETHYL  KETONE 

1.22 

Mh'l'HYL  ISOBUTYL  KETONE 

0.56 

MaTrmM,T7MF 

0.18 

PERCHDDRDETHYLENE 

0.70 

PHENOL 

0.44 

PRDPIONAIDEHYDE 

0.42 

PROPIONIC  ACID 

0.40 

PROPYLENE  OXIDE 

0.50 

STYRENE 

1.20 

TOLUENE 

2.04 

TOLUENE  DIISOCYANATE 

0.66 

TRICHLDRDBENZENES 

0.01 

TRICHDDRDETHYIENE 

0.43 

VINYL  CHIDRIDE 

0.42 

OTHER 

76.67 

Subtotal 

99.96 
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SIC  2919  Other  Primary  Steel  IncJustries 
IM  Contaminant    Per  Cent 


ARSENIC 

'■:i:      0.08 

CAEMIUM 

■    .      ■      0.09 

CHRCMIUM 

0.76 

FIUORIDE 

0.07 

MANGANESE 

5.06 

NICKEL 

0.20 

LEAD 

1.63 

SKTFTmM 

0.01 

SHIGA 

0.85 

OTHER 

91.25 

Siobtotal 

100.00 

VDC  Contaminant  Per  Cent 

1,3-BUTADIENE  0.18 

ACRYIJDNITRILE  0.19 

BENZENE  2.16 

BUTYL  ACRYIATE  0.04 

BUTYL  CELLDSOLVE  2.09 

CARBON  TETRACHLDRIDE  0.06 

CHLDRDFORM  0.01 

ETHYIEENZENE  0.30 

ETHYLENE  l-^S 

ETHYLENE  DIGHDDRIDE  0.50 

ETHYLENE  OXIDE  0.01 

PORMAIDEHYDE  0.54 

METHYL  ETHYL  KETONE  2.56 

METHYL  ISOBUTYL  KETONE  1.79 

NAPIHALENE  0.03 

PHENOL  0.33 

STYRENE  0 .  89 

TOUJENE  20.36 

TRICHLDRDBENZENES  0 . 01 

VINYL  CHDDRIDE  0.32 

OTHER  65.88 

Subtotal  100.01 
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POUNCRIES 


SIC  2941  Iron  Fcundries 


FM  Contaminant  ;,-;,  Per  Cent 


ARSENIC 

0.01 

CAEKCUM 

0.01 

CHROOUM 

0.08 

MANGANESE 

3.24 

NICKEL 

0.10 

LEAD 

0.38 

SILICA 

13.58 

OTHER 

82.59 

Subtotal 

99.99 

VOC  Contaminant 

Per  Cent 

1,1, l-TRICHDDROETHANE 

0.05 

1,3-BOTADIENE 

0.08 

ACRYIONl'i'RILE 

0.07 

BENZENE 

1.42 

BUTYL  ACRYIATE 

0.05 

BUTYL  CELLDSOLVE 

4.04 

CARBON  'i'E'i'KACHLDRIDE 

0.08 

rvTjr\^m\7P. 

0.04 

CKI,U)SOLVE  ACETATE 

0.04 

CHIDRDFORM 

0.06 

EPICHIDRDHYDRIN 

0.05 

ETHYLBENZENE 

0.41 

ETHYLENE 

0.36 

ETHYLENE  DIBRCMIDE 

0.05 

ETHYLENE  DICHIDRIDE 

0.07 

ETHYLENE  OXIDE 

0.05 

FORMALDEHYDE 

0.21 

METHYL  ETHYL  KETONE 

0.48 

METHYL  ISOBUTYL  KETONE 

0.29 

NAPTHATENE 

0.02 

PERCHLDROETHYLENE 

0.09 

PHENOL 

3.13 

PROPIONAIDEHYDE 

0.05 

PROPIONIC  ACID 

0.05 

PROPYLENE  OXIDE 

0.06 

STYRENE 

0.15 

TOILIENE 

23.97 

TOLUENE  DIISOCYANATE 

0.08 

TRICHLDROETHYLENE 

0.05 

VINYL  CHLDRIDE 

0.05 

OTHER 

64.34 

Subtotal 

99.94 
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SIC  2971  C3opper  and  Copper  Alloy 
Rolling,  Casting  and  Extruding 


IM  Contaminant 

•  Per  Cent 

ARSENIC 

••    -     ••     0.03 

CAEMIUM 

0.01 

OiRCMrUM 

0.02 

MANGANESE 

0.71 

NICKEL 

0.02 

LEAD 

0.37 

Sn,TCA 

3.60 

TIN 

0.01 

OTHER 

95.25 

Subtotal 

100.02 

VDC  Contaminant 

Per  Cent 

1,1, 1-TRICHLORDEIHANE 

6.93 

BENZENE 

0,31 

BUTYL  CKUnSOLVE 

3.33 

ETHYLBENZENE 

0.28 

MhTH^L  EiHI^L  Kh'iONE 

0.62 

ML'i'HYL  ISOBOTYL  KETONE 

0.18 

NAPTTTATKNE 

0.01 

i^l::RCHU)WDEniYIENE 

18.30 

TOIUENE 

22.03 

TRICHLDRDETHYLENE 

6.55 

OTHER 

41.45 

Subtotal  1 

99.99 

SIC  2999  Other  Rolled,  Cast  and  Extnided 
Ncn-Fercous  Metal  Products 


IM  Contaminant 

Per  Cej 

ARSENIC 

0.06 

CAEMIUM 

0.05 

CHRCMIUM 

1.85 

MANGANESE 

0.61 

NICKEL 

0.53 

LEAD 

0.95 

SILICA 

4.80 

TIN 

0.06 

OTHER 

91.08 

Subtotal 

99.99 

VDC  Contaminant 

PERCHLDRDETHYLENE 
Subtotal 


Per  Cent 

100.00 
100.00 
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NCN-irlKWLXJS  MEIAL  SMELnNG  AND  REFINING 

SIC  2959     Other  Primary  Smelting  and 
Refining  of  Ncn-Ferrous  Metal 


IM  Contaminant    - 

Per  Cent 

ARSENIC 

25.48 

CAEMIUM 

1.46 

CHRC3-1IUM 

0.22 

FITORIDE 

1.47 

MERCURY 

0.04 

MANGANESE 

0.32 

NICKEL 

0.09 

LEAD 

12.25 

SFTENIUM 

0.23 

snicA 

2.25 

TIN 

0.32 

OTHER 

55.86 

Subtotal 

99.99 

VOC  Contaminant 

BENZENE 

ETHYLENE 

FORMALDEHYDE 

TOIIJENE 

OTHER 

Subtotal 


Per  Cent 

4.01 

0.01 

7.98 

2.00 
85.99 
99.99 


mnOR  VEHICLE  INDCSTRIES 

SIC  3231    Ifctor  Vehicles  and  Motor 
V^iicle  Parts  and  Accessaries  Surface 
Coating 


H-I  Contaminant 

Per  Cent 

CAEMILT^ 

0.05 

CHRCMUM 

0.05 

NICKEL 

0.05 

OTHER 

99.85 

Subtotal 

100.00 

VDC  Contaminant 

BENZENE 

ETHYIBENZENE 

TOIIJENE 

OTHER 

Subtotal 


Per  Cent 

1.04 

8.68 

20.12 

70.16 

100.00 
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SIC  3254  Motor  Vehicle  Steering  and 
Suspensicai  Parts 


FM  Ccntaminant 


Per  Cent 


VDC  Contaminant 


OTHER 
Subtotal 


99.99 
99.99 


1,1,1-TRICHIOROETHANE  1.89 

1,3-BOTADIENE  0.96 

ACRYIDNITRILE  0 .  01 

BENZENE  0.39 

BOTYL  CELDDSOLVE  0.76 

ETHYLBENZENE  2.33 

ETHYLENE  2.29 

ETHYLENE  DICHIDRIDE  0.02 

FORMAUDEHYDE  0.44 

METHYL  ETHYL  KETONE  0.17 

METHYL  ISOBUTYL  KETONE  0.04 

PERCHLDROETHYLENE  0 .  63 

HffiNOL  0.01 

PRDPIONAUDEHYDE  0.01 

STYRENE  0.04 

TOUJENE  12.76 

TRICHDDRDETHYLENE  1.79 

VINYL  CHLORIDE  0.01 

OTHER  75.45 

Subtotal  100.00 


XNQRGANIC  CHEMICaL  MANUEACTORING 


SIC  3511  Domestic  day  Products 
IM  Contaminant    Per  Cent 


VOC  Contaminant 


Per  Cent 


ARSENIC 

0.01 

CATKUM 

0.01 

CHRCMIUM 

0.17 

FIUORIDE 

0.03 

MANGANESE 

0.07 

NICKEL 

0.06 

LEAD 

0.17 

SKTFNIUM 

0.06 

SIUCA 

7.67 

ct:her 

91.76 

Subtotal  1 

100.01 

BENZENE 

ETHYLBENZENE 

ETHYLENE 

PORMALCEHYDE 

NAPIHALENE 

TOLUENE 

OHiER 

Subtotal 


3.76 
3.29 
0.08 
2.96 
0.10 
5.59 
84.21 
99.99 
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SIC  3512     Inported  day  Products 

FM  Contaminant    Per  Cent 


ARSENIC 

.:.•!.■  0.01 

CAEMIUM 

'  -  ■■■   0.01 

CHRCMIUM 

0.12 

FLDORTDE 

0.02 

MANGANESE 

0.04 

NICKEL 

0.04 

lEAD 

0.11 

SFTFNIUM 

0.04 

Sn.TCA 

6.72 

OTHER 

92.90 

Subtotal 

100.01 

VDC  Contaminant 

Per  Cent 

1,3-BUrADIENE 

0.01 

ACRYTONI'I'RTTK 

0.01 

Bmzwb: 

0,09 

BUIYL  (TTinSOLVE 

2.63 

ETHYTRKNZENE 

0.81 

ETHYLENE 

0.08 

ETHYLENE  DIERCMDE 

0,09 

ETHYLENE  DICHLDRIDE 

0.02 

FORMAUDEHYDE 

0.07 

ML'iHYL  EIHYL  KETONE 

2.56 

ML"i'H!i:L  iSOBLTi^L  Kh'iONE 

1.62 

FHENOL 

0.01 

STYRENE 

0.04 

TOIUENE 

25.71 

VINYL  CHLORIDE 

0.01 

OTHER 

66.22 

Subtotal 

99.98 

E  -  25 


SIC  3521    HycJraulic  Oeanent 

FM  Contaminant    Per  Cent 


CAEMIUM 

•  ':,:{.   0.01 

CHRCMIUM 

-  ■■  .  ■  0.22 

FTUORIDE 

0.12 

MANGANESE 

0.02 

NICKEL 

0.03 

LEAD 

0.83 

SETramJM 

0.01 

SILICA 

5.48 

OTHER 

93.29 

Subtotal 

100.01 

VDC  Contaminant 

Per  Cent 

1,1,  l-TRICHLDRDETHANE 

0.25 

1,3-BOTADIENE 

0.39 

ACRYTDNl'l'Rnf; 

0.34 

BENZENE 

2.60 

BUTYL  ACRYIATE 

0.25 

TTTTVT,  rvrjnfim\7P. 

0.26 

CARBON  I'E'i'RACHLDRIDE 

0.38 

cFT,Ta=;oTVF: 

0.21 

CKIJOSULVE  ACL'IME 

0.21 

CHLOROFORM 

0.30 

EPICHLDROHYDRIN 

0.25 

eihyt,r?:nzene 

0.56 

ETHYLENE 

2.15 

ETHYLENE  DIBRCMIDE 

0.25 

ETHYLENE  DICHIORIDE 

0.34 

ETHYLENE  OXIDE 

0.22 

FORMAIDEHYDE 

4.95 

MLTHVL  hTHVL  KL'iONE 

0.72 

ML'iHYL  ISOBUIYL  KETONE 

0.33 

NAFTHATKNE 

0.27 

PKRCHUDROEIHYLENE 

0.41 

IHENOL 

0.26 

PROPIONAIDEHYDE 

0.25 

PROPIONIC  ACID 

0.24 

PROPYLENE  OXIDE 

0.29 

STYRENE 

0.71 

TOLUENE 

1.76 

TOLUENE  DIISOCYANATE 

0.39 

TRICHLOROBENZENES 

0.01 

TRICHLOROETHYIENE 

0.25 

VINYL  CHLORIDE 

0.25 

OTHER 

79.91 

Subtotal 

99.96 

3571    Abrasives 


E  -  26 


EM  Contaminant 

Per  Cent 

ARSENIC 

;.-,',  0.47 

CAEMIUM 

0.15' 

CHRCMIUM 

0.05 

EUJOPTnE 

0.10 

MERCURY 

0.02 

MANGANESE 

0.14 

NICKEL 

0.05 

LEAD 

1.40 

STJFKHM 

0.01 

STTiTGA 

3.47 

TIN 

0.06 

OTHER 

94.08 

Subtotal! 

100.00 

VDC  Contaminant 

Per  Cent 

1,1/  1-IRICHIORDEIHANE 

0.02 

1,3-BOTADIENE 

0.13 

ACRYTnNl'I'KTTF: 

0.13 

BENZKNE 

1.26 

BUTYL  ACRYIATE 

0.05 

BUTYL  CETJnSOLVE 

0.02 

CARBON  'i'LTRACHLDRIDE 

0.07 

(TTJnSOLVE 

0.02 

rVJJDf^TTJF.   ACFTfATR 

0.02 

CHDDRDPORM 

0.03 

EPICHIORDHYnRIN 

0.02 

ETHYTRENZENE 

0.84 

ETHYLENE 

1.09 

ETHYLENE  DIHHCMIDE 

0.74 

ETHYLENE  DICHIDRIDE 

0.30 

ETHYLENE  OXIDE 

0.03 

FORMALDEHYDE 

0.37 

ML'i'HYL  ETHYL  KETONE 

3.47 

ML'iHYL  ISOBUTYL  KETONE 

0.06 

NAFTHATENE 

0.02 

PKkOlLUROEIHYIENE 

0.04 

rHENOL 

0.20 

ERDPIONAIDEHYDE 

0.02 

PROPIONIC  ACID 

0.02 

PROPYLENE  OXIDE 

0.03 

STYRENE 

0.54 

TOLUENE 

12.36 

TOLUENE  DIISOCYANATE 

0.03 

TRICHIOROBENZENES 

0.01 

TRICHIDROETHYLENE 

0.02 

VINYL  CHLORIDE 

0.19 

OTHER 

77.84 

STJbtotal 

99.99 

E  -  27 


SIC  3581 

T.iinft 

m.  Contaminant 

Per  Cent 

ARSENIC 

> 

■:,1>  0.01 

CAEMIUM 

.  ; 

-  ■  0.01' 

CHRCMIUM 

0.17 

FLUORIDE 

0.03 

MANGANESE 

0.07 

NICKEL 

0.06 

LEAD 

0.17 

sEiramM 

0.06 

SHIGA 

7.67 

OTHER 

91.76 

Subtotiil 

100.01 

SIC  3711 

Industrial  Inorya 

(NEC) 

FM  Contaminant 

Per  Cent 

CAHGUM 

0.04 

CHRCMIUM 

0.01 

MANGANESE 

0.03 

NICKEL 

0.04 

LEAD 

0.02 

SFTFTTIUM 

0.01 

SILICA 

1.62 

TIN 

0.01 

OTHER 

98.23 

Subtotal 

100.01 

VDC  Contaminant 

BENZENE 

EIHYLBENZENE 

ETHYLENE 

PORMAIDEmrDE 

TOLUENE 

OTHER 

Sxobtotal 


Per  Cent 

1.38 

0.81 

1.67 
15.19 

1,48 
79.46 
99.99 


VDC  Contaminant 

Per  Cent 

1,1, 1-TRICHIORDETHANE 

0.54 

1,3-BUTADIENE 

0.56 

ACRYIONl'i'RIIE 

0.23 

BENZENE 

3.19 

BUTYL  ACRYIATE 

0.06 

BUTYL  nrrJOSOLVE 

0.01 

CARBON  Th'i'KACHDDRIDE 

0.08 

CHLDRDFORM 

0.02 

EPICHLDROHYDRIN 

0.01 

EIHYTRT\NZENE 

0.72 

ETHYLENE 

4.21 

ETHYLENE  DIBRCMIDE 

0.01 

ETHYLENE  DICHDDRIDE 

2.10 

ETHYLENE  OXIDE 

0.02 

FORMALDEHYDE 

0.11 

ML'iHYL  EiHYL  KETTONE 

0.49 

METHYL  ISOBOTYL  KETONE 

0.07 

PEKCHIOROETHYLENE 

0.19 

PHENOL 

0.37 

PROPIONAIDEHYDE 

0.01 

PROPIONIC  ACID 

0.01 

PROPYLENE  OXIDE 

0.01 

STYRENE 

5.74 

TOLUENE 

0.85 

TOLUENE  DIISOCYANATE 

0.01 

TRICHLDROErHYLENE 

0.51 

VINYL  CHLORIDE 

0.38 

OIHER 

79.52 

Subtotal  1 

100.03 

E  -  28 

SIC  3721  Cheanical  Fertilizer  and 
Fertilizer  Materials 


IM  Contaminant 

Per  Cent 

ARSENIC 

'.'"  0.05- 

CAEMIUM 

0.05 

CHRCMIUM 

0.03 

FIUORIDE 

0.01 

MANGANESE 

0.09 

NICKEL 

0.04 

lEAD 

0.16 

SKTKNIUM 

0.03 

STT.TCA 

4.04 

TIN 

0.02 

OIHER 

95.46 

Subtotal 

99.98 

VDC  Contaminant 

Per  Cer 

1,1, l^IRICHIDRDETHANE 

0.01 

1,3 -BUTADIENE 

0.02 

ACRYlJDNiTKnK 

36.80 

BENZENE 

2.73 

BUTYL  ACRYIATE 

0.07 

BUTYL  rKTiTOSOLVE 

0.01 

CARBON  TEi'RACHDDRIDE 

0.10 

rPT.TnsoTVR 

0.01 

CK]  ,1  r«OLVE  ACh'iATE 

0.01 

CHDDRDPORM 

0.03 

EPiCHIDRDHYDRIN 

0.01 

ETHYTHENZENE 

0.19 

ETHYLENE 

2.31 

ETHYLENE  DlHtOttDE 

0.01 

ETHYLENE  DICHDDRIDE 

0.69 

ETHYLENE  OXIDE 

0.03 

FORMALDEHYDE 

0.78 

MhTH!i!L  EiHVL  KL'iONE 

0.55 

ML'i'HYL  ISOBUTYL  KETONE 

0.08 

MAPTTOTFTTO 

0.01 

jhekuujjwdeihylene 

0.02 

HIENOL 

0.45 

ERDPIONALDEHYDE 

0.01 

PROPIONIC  ACID 

0.01 

PROPYLENE  OXIDE 

0.02 

STYRENE 

25.44 

TOLUENE 

0.73 

TOrUENE  DIISOCYANAIE 

0.02 

TRICHLDRDETHYLENE 

0.01 

VINYL  CHLORIDE 

0.44 

XYLENE  BASE  ACIDS 

0.01 

OIHER 

28.35 

Subtotal 

99.96 

E  -  29 


PETODIECM  RKJ^'INEHTKS 


SIC  3611  Petnaleum  Refineries 


IM  Contaminant  ;,-^'Per  Cent 


ARSENIC 

0.01 

CHRCMIUM 

0.14 

FLUORIDE 

0.02 

MANGANESE 

0.03 

NICKEL 

1.51 

LEAD 

0.07 

SKTENIUM 

0.13 

STTJCA 

8.57 

OTHER 

89.51 

Subtotal 

99.99 

VDC  Contaminant 

BENZENE 

ETHYIBENZENE 

ETHYLENE 

PORMAIDEHYDE 

NAPIHALENE 

STYRENE 

TOUJENE 

OTHER 

Subtotal 


Per  Cent 

2.60 
0.36 
0.01 
6.69 
0.07 
0.01 
1.75 
88.51 
100.00 


ASTHAIir  PAVING  INDUSTRY 

SIC  3699  Asphalt  Paving  Industry 

IM  Contaminant    Per  Cent 


ARSENIC 

0.01 

CAEMIUM 

0.01 

CHRCMIUM 

0.15 

FLUORIDE 

0.02 

MANGANESE 

0.07 

NICKEL 

0.05 

LEAD 

0.15 

SKTFNIUM 

0.06 

SILICA 

6.73 

OTHER 

92.75 

Subtotal 

100.00 

VDC  Contaminant 

Per  Cent 

1,3-BOTADIENE 

0.01 

ACRYIONl'i'RILE 

0.01 

BENZENE 

2.95 

CARBON  TEl'RACHLDRIDE 

0.01 

CHLOROFORM 

0.01 

EIHYT, BENZENE 

1.25 

ETHYLENE 

0.03 

ETHYLENE  DICHIDRIDE 

0.01 

PORMAIDEHYDE 

5.53 

Mh'iHYL  ETHYL  KETONE 

0.01 

ML'IHYL  ISOBUTYL  KETONE 

0.01 

MAPTHATFMF 

1.02 

PKHCHLOROEIHYIENE 

0.01 

PROPYLENE  OXIDE 

0.01 

STYRENE 

0.01 

TOLUENE 

4.28 

TOLUENE  DTISOCYANATE 

0.01 

OTHER 

84.81 

Subtotal 

99.98 

E  -  30 


OE^GANIC  CHEJIIC2UL  ERDDOdS 

SIC  3712  Industrial  Organic  Qiemicals 
(NEC) 


FM  Contaminant- 

•  Per  Ceri 

ARSENIC 

0.01 

CAEMIUM 

0.04 

CHPCMIUM 

0.02 

MERCURY 

0.16 

MANGANESE 

0.09 

NICKEL 

0.06 

LEAD 

0.27 

SFTKNIUM 

0.02 

STT.TCA 

4.45 

TIN 

0.05 

OTHER 

94.84 

Subtot;^! 

100.01 

VOC  contaminant 

Per  Cent 

1, 1,  l-TRICHKJROEIHANE 

0.01 

1,3 -BUTADIENE 

2.08 

acryidni'ikttf: 

4.11 

BENZENE 

4.32 

BUTYL  ACRYIATE 

0.08 

BUTYL  CKTiTDSOLVE 

0.01 

CARBON  TEIHACHDDRIDE 

0.12 

CHDDRDFORM 

0.02 

EPiCHTJDRDHYDRIN 

0.01 

ETHYTRFNZENE 

0.50 

ETHYLENE 

16.89 

ETHYLENE  DIHHCMIDE 

0.01 

ETHYLENE  DICHIORIDE 

2.36 

ETHYLENE  OXIDE 

0.03 

PORMAIDEHYDE 

0.05 

MEiHYL  EIHYL  KETONE 

2.72 

ML'iH!i!L  ISUBUiYL  Kh'iONE 

0.09 

PERCHDDRDETHYLENE 

0.01 

HffiNOL 

0.87 

PROPIONAIDEHYDE 

0.01 

PROPIONIC  ACID 

0.01 

PROPYLENE  OXIDE 

0.02 

STYRENE 

5.50 

TOIDENE 

0.80 

TOIUENE  DIISOCYANATE 

0.01 

TRICHLDRDETHYLENE 

0.01 

VINYL  CHDDRIDE 

3.16 

XYIENE  BASE  AC1I« 

0.01 

OTHER 

56.17 

Subto+val 

100.00 

SIC  3731  Plastics  and  Synthetics  Resins 

IM  Contaminant    Per  Cent 


CAEMIUM 

0.04 

CHROMIUM 

0.02 

MANGANESE 

0.12 

NICKEL 

0.03 

LEAD 

0.01 

SFTENIUM 

0.02 

SILICA 

5,95 

TIN 

0.03 

OTHER 

93.80 

Sv±tot;^l 

100.02 

VOC  Contaminant 

BENZENE 

PORMAIDEHYDE 

TOLUENE 

OTHER 

Subtotal 


Per  Cent 

0.70 

1.41 

0.35 

97.55 

100,01 


E  -  31 


SIC  3751  Paint  and  Vaznish 
IM  Contaminant    "Per  Cent 


CAEMIUM 

•  ••    -         0.03 

CHRCMIUM 

0.02 

MANGANESE 

0.12 

NICKEL 

0.02 

LEAD 

0.01 

SFrn^IUM 

0.01 

SIUCA 

5.93 

TIN 

0.03 

OTHER 

93.81 

Subtotal 

99.98 

SIC  3791  Printing  Ink 

IM  Contaminant     Per  Cent 


CAEMIUM 

0.04 

CHRCMIUM 

0.02 

MANGANESE 

0.12 

NICKEL 

0.03 

LEAD 

0.01 

SKTFOTUM 

0.02 

STT.TCA 

5.95 

TIN 

0.03 

OTHER 

93.80 

Subtotj=il 

100.02 

VDC  Contaminant 

Per  Cent 

1,3-BOTADIENE 

10.45 

EENZHffi 

0.01 

ETHYLENE 

0.01 

FORMALDEHYDE 

0.01 

Mh'lHyL  K1H!s:L  KL'iONE 

7.61 

PHENOL 

5.36 

STYRENE 

6.97 

TOLUENE 

31.70 

OTHER 

37.84 

Subtotal 

99.96 

VOC  Contaminant 

Per  Cent 

ETHYLBENZENE 

3.00 

Mh'i'HVL  EiHYL  Kh'lONE 

5.00 

OIHER 

92.00 

Svibtoti^l 

100.00 

E  -  32 


SIC  3792  Adhesives 


EM  Contaminant 

Per  Cer 

CAEMIUM 

■  ■'   0.05 

CHRCMIUM 

V''  0.04- 

MANGANESE 

0.04 

NICKEL 

0.04 

SETRTIUM 

0.01 

SIUCA 

1.98 

TIN 

0.01 

OIHER 

97.84 

Sx±)total 

100.01 

vex:  Contaminant 

Per  Cent 

l,3-Bm:ADIENE 

0.07 

ACRYIDNllMLE 

0.07 

BENZENE 

0.82 

BOTYL  ACRYLATE 

0.02 

BOTYL  rKTTOSOLVE 

0.09 

CARBON  'i'h'l'KACHIDRIDE 

0.02 

CELLDSOLVE 

0.57 

rrrjr^qnTVF  acktatr 

0.57 

ETHYLBENZENE 

0.12 

ETHYIENE 

0.69 

ETHYLENE  DIBRCMIDE 

0.01 

ETHYLENE  DIOfflDRIDE 

0.19 

ETHYLENE  OXIDE 

0.57 

FORMALDEHYDE 

0.21 

ML'i'HYL  ETHYL  KETONE 

1.08 

Mh'i'HYL  ISOBUTYL  KETONE 

0.15 

NAPTHATraJE 

0.01 

EHENOL 

9.03 

STYRENE 

0.35 

TOUJENE 

1.32 

VINYL  CHDDRIDE 

0.12 

OTHER 

83.89 

Subtotal 

99.97 

SIC  3799  Other  Qiemical  Products  (NEC) 

FM  Contaminarit     Per  Cent 


CAEMIUM 

0.03 

CHRCMIUM 

0.05 

MANGANESE 

0.05 

NICKEL 

0.07 

LEAD 

0.01 

SEr,KNIUM 

0.05 

SILICA 

3.05 

TH^I 

0.01 

OTHER 

96.70 

Subtotal 

100.02 

VOC  Contaminant 

Per  Cent 

ACRYLDNITRILE 

0.11 

BENZENE 

0.74 

BUTYL  ACRYIATE 

0.01 

CARBON  I'h'i'RACHLDRIDE 

0.02 

ETHYT, BENZENE 

0.04 

ETHYLENE 

0.91 

ETHYLENE  DICHIORIDE 

0.17 

ETHYLENE  OXIDE 

0.08 

FORMALDEHYDE 

0.03 

ME'i'HYL  ETHYL  KETONE 

0.13 

ML'iHYL  ISOBUTYL  KETONE 

0.02 

FHENOL 

10.44 

PROPYLENE  OXIDE 

0.02 

STYRENE 

0.32 

TOLUENE 

0.19 

VINYL  CHIDRIDE 

0.12 

XYLENE  BASE  ACIDS 

0.02 

OTHER 

86.62 

SxJbtotal 

99.99 

E  -  33 


EIECmrC  POWER  SYSTEMS 

SIC  4911  Electric  BDWsr  Systons  -  C3oal 
and  Heavy  Oil  FxxpeQ.  Generating  Stations 


FM  Contaminant 

Per  Cent 

ARSENIC 

0.05 

CHRCMIUM 

0.06 

MERCURY 

0.01 

MANGANESE 

0.04 

NICKEL 

0.25 

LEAD 

0.17 

SFTJ\NIUM 

0.02 

STT.TCA 

20.06 

TIN 

0.01 

OTHER 

79.32 

Subtotal 

99.99 

VDC  Contaminant  Per  Cent 

1,3-BUTADIENE  0.28 

BENZENE  0.37 

ETHYLBENZENE  5.88 

ETHYLENE  1.87 

FORMALDEHYDE  14.19 

PRDPIONAIDEHYDE  0 . 01 

TOUJENE  2.96 

OIHER  74.39 

Subtotal  99.95 


WASTE  DISPOSAL  INDUSTRIES 

SIC  4999  Hospital  Incinerators 

EM  Contaminant    Per  Cent 


ARSENIC 

0.01 

CAEMIUM 

0.06 

CHRC^aUM 

0.02 

MERCURY 

2.67 

MANGANESE 

0.04 

NICKEL 

0.01 

LEAD 

4.12 

TIN 

0.44 

OTHER 

92.64 

Subtotal         100.01 


E  -  34 


Incinerators 

m.  Contaitiinant 

'Per  Cent 

APSENIC 

'- '      Oc02' 

CAEMIUM 

0.09 

CHRCMIUM 

0,03 

0.02 

MERCURY 

3.99 

MANGANESE 

0.07 

NICKEL 

0.01 

LEAD 

6.15 

STTJCA 

0.05 

TIN 

0.66 

OTHER 

88.93 

Subtotal 

100.02 

VDC  Contaminant 

BENZENE 
ETHYLENE 
COHER 
Subtotal 


SIC  4999  Sewage  SlucJge  Incinerators 

IM  Contaitiinant    Per  Cent 


CAEMIUM 

0.83 

CHRCMIUM 

0.58 

MANGANESE 

0.21 

IZAD 

2.94 

SKTITNIUM 

0.15 

snzcA 

4.77 

TIN 

1.39 

OTHER 

89.13 

Subtotal 

100.00 

ERY  aFT^NINS  ESTABI 

J[SHME^^ 

VDC  Contaminant 

BENZENE 
ETHYLENE 
OTHER 
Subtotal 


SIC  9721  Dry  Cleaners 

IM  Contaminant    Per  Cent 


ARSENIC 

1.18 

CAEt-nUM 

0.38 

CHRCMIUM 

0.13 

FIDOPTnP; 

0.24 

MERCURY 

0.05 

MANGANESE 

0.34 

NICKEL 

0.13 

LEAD 

3.50 

SKTFNILH 

0.03 

Sn.TCA 

8.67 

TIN 

0.16 

OTHER 

85.19 

Subtotal 

100.00 

VDC  Contaminant 

BENZENE 

ETHYLBENZENE 

FORMAIDEHYDE 

NAPIHAIENE 

PERCHIDROETHYLENE 

PHENOL 

TOIDENE 

COHER 

Subtotal 


Per  Cent 

7.28 

9.27 
83.43 
99.98 


Per  Cent 

7.70 

8.70 

83.60 

100.00 


Per  Cent 

0.03 

0.16 

0.07 

0.16 
13.23 

4.93 

0.25 

81.20 

100.03 
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F.l   INTRODUCTION 


This  appendix  (describes  the  dispersion  modelling  work  performed 
to  estimate  gi'ids  for  annual  average  concentration  of  air 
emissions  in 'various  parts  of  Ontario. 

F.1.1   Organization  of  the  Work 

The  principal  steps  in  the  work  are: 
.  selection  of  a  model; 

.  selection  of  locations  for  which  the  model  is  run; 
.  compilation  of  required  meteorological  data  for  each 

location; 
.  specification  of  source  characteristics;   and 
.  quality  assurance  checks  of  the  model  results. 

The  approach  taken  to  each  of  these  steps  is  outlined  in  the 
following  sections. 

F.1.2   Acknowl edg ement 
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L.  Bradizza. 

We  also  appreciate  the  assistance  of  staff  of  the  Canadian 
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F.2   SELECTION  OF  A  DISPERSION  MODEL 


F.2.1  Proposed  Dispersion  Model 

Our  original  proposal  described  a  regulatory  dispersion  model 
presently  in  use  in  the  Province  of  Alberta.   This  model  was  a 
variant  of  the  existing  Ontario  regulatory  model  (a  Pasguill  and 
Gifford  type) .   The  multi-source,  receptor  version  of  the  Alberta 
model  was  considered  superior  in  many  respects  to  the  existing 
Ontario  regulatory  model  and  thus  to  be  preferred. 

Due  to  demand  generated  by  widespread  use  of  Pasquill-Gifford 
models,  Environment  Canada  has  developed  the  STAR  (STability 
ARray)  program  to  generate  the  simplified  climatological  input 
data  used  by  these  models.   Typical  input  consists  of  frequency 
distributions  of  wind  direction,  wind  speeds,  and  the  Pasquill 
stability  classes  (A,  B,  C,  etc.). 

F.2. 2   Ministry  of  the  Environment's  Models 

In  discussions  with  the  Ontario  Ministry  of  the  Environment's 
Atmospheric  Model  Development  Unit,  during  the  early  stages  of 
this  project,  it  was  agreed  that  use  of  models  currently  under 
development  in  Ontario  would  be  preferred.   The  models  under 
development  were  far  superior  to  the  Pasquill-type  models.   The 
only  relative  shortcomings  of  the  newer  models  were  in  the  ease 
of  use.   For  example,  meteorological  input  data  are  in  the  form 
of  detailed  hourly  meteorological  observations.   This  represents 
a  quantum  leap  in  the  quantity  of  data  to  be  handled  and  gave 
rise  to  considerable  complications  in  data  transfer  and  data 
management . 

A  development  version  of  the  Ministry's  multi-source,  (GAS)  model 
was  released  to  the  study  team  in  early  June.   This  helped 
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facilitate  data  conversion  program  development.  A  final 

"pre-public-release"  version  was  released  to  the  study  team  on 

June  21,  1988.'  That  version  was  used  to  execute  the  modelling 
runs.      •  ''^',. 

Some  minor  modification  to  the  meteorological  parameter 
calculation  program  (MET)  provided  by  the  Ministry  was  undertaken 
to  flag  additional  meteorological  cases  as  invalid  (so  these 
cases  would  not  be  calculated) .   Program  and  data  conversion 
development  are  described  in  more  detail  in  subsequent  sections. 

A  full  discussion  of  the  new  models  is  presented  in  Appendix  H  of 
the  Clean  Air  Program  Discussion  Paper  issued  by  the  Ontario 
Ministry  of  the  Environment  in  late  1987  (Environment  Ontario, 
1987)  . 


F.3   SELECTION  OF  LOCATIONS 


F . 3 . 1   Introduction 

To  determine  the  public  benefits  of  the  proposed  revisions  to 
Regulation  308,  it  is  necessary  to  determine  the  dispersion  of 
emissions  of  each  of  the  96  specified  contaminants  emitted  by 
each  of  the  over  3,500  sources  identified.   In  principle,  each 
source  must  be  modelled  separately  because  the  source 
characteristics  and  meteorological  conditions  are  unique.   In 
practice  that  is  impossible  due  to  lack  of  data  and  the  magnitude 
of  the  effort  it  would  require. 

The  availability  of  climatological  information  allows  only  a 
limited  number  of  sites  to  be  modelled;   those  for  which  the 
required  input  data  can  be  compiled.   The  dispersion  patterns 
modelled  for  those  locations  are  then  applied  to  all  sources 
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within  the  area  to  which  those  cliitiatological  characteristics 
apply.   This,  in  turn,  requires  that  standard  characteristics  be 
assumed  for  the  emission  sources.   The  use  of  standard  source 
characteristics  is  not  a  major  concern  since  the  required 
information  is,  with  very  rare  exceptions,  not  available  in  any 
case.   The  source  characteristic  assumptions  are  discussed  in 
Section  F.5  below.   Here  we  review  the  procedures  used  to 
identify  the  locations  with  different  climatological 
characteristics . 

F.3.2   Climatological  Data  Requirements 

The  Ministry's  GAS  model  requires  input  climatological  data  from 
both  hourly  observing  stations  and  upper  air  stations.   There  are 
no  locations  within  the  study  area  for  which  both  hourly 
observations  and  upper  air  data  are  available.   To  get  the  model 
input  required  upper  data  from  U.S.  and  Quebec  upper  air  stations 
must  be  matched  with  data  from  Ontario  hourly  observing  stations. 

The  process  of  selecting  locations  to  be  modelled  starts  with  the 
list  of  hourly  observing  stations  in  Ontario  and  eliminates  those 
that  have  similar  climatological  conditions.   Appropriate  upper 
air  data  are  then  matched  to  these  remaining  locations. 

The  Ministry's  MET  and  GAS  models  require  hourly  air  temperature, 
cloud  opacity,  wind  speed  and  wind  direction  observations  in 
addition  to  a  daily  snow  on  the  ground  indicator  and  mixing 
height,  which  is  derived  from  upper  air  temperature  profiles. 
Thus  the  starting  point  for  the  selection  of  locations  is  the  32 
hourly  observing  stations  in  Ontario  in  the  climatological 
archive  of  Environment  Canada,  Atmospheric  Environment  Service. 
Most  of  these  stations  have  a  data  record  extending  for 
approximately  3  0  years.   This  is  the  typical  period  of  record 
considered  necessary  for  a  useful  climatological  data  set. 
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One  year  is  composed  of  8760  hours,  so  the  analysis  of  four 
hourly  weather  elements  for  a  30  year  period  for  all  Ontario 
stations  would'  require  analysis  of  over  3  0  million  observations. 
This  was  well' 'beyond  the  scope  of  the  project. 

F.3.3   Approach  to  the  Selection  of  Locations 

To  reduce  the  number  of  locations,  and  hence  the  computational 
work  load,  in  a  scientifically  credible  manner  it  is  necessary 
to  develop  criteria  to  determine  when  stations  have  similar 
climatological  characteristics  from  the  perspective  of  the 
dispersion  model. 

An  examination  of  the  Ontario  Ministry  of  the  Environment's  model 
indicated  that  the  development  of  dispersion  contours  around  a 
source  depends  heavily  on  the  wind  speed  and  direction  of  the 
source  area.   Accordingly,  we  began  by  comparing  mean  annual 
statistics  of  wind  speed  and  direction  for  the  32  available 
stations. 

Several  stations  are  located  very  close  to  one  another,  and  so 
allow  us  to  define  criteria  for  similar  climatological 
characteristics.   For  example,  Toronto  Downsview  Airport 
(latitude  43°  45*  N,  longitude  79°  29*  W)  or  Toronto 
International  Airport  (latitude  43°  40 '  N,  longitude  79°  38'  W) 
could  be  used  to  represent  conditions  in  the  Toronto  area. 
Similarly,  either  Ottawa  Rockcliffe  Airport  (latitude  45°  27*, 
longitude  75°  38'  W)  or  Ottawa  International  Airport  (latitude 
45°  19   N,  longitude  75°  40 '  W)  could  be  used  to  represent 
conditions  in  the  Ottawa  area.   It  follows  that  if  other 
geographically  proximate  stations  had  wind  statistics  as  similar 
to  one  another  as  these  station  pairs,  then  they  too  could  be 
considered  interchangeable  for  the  purposes  of  this  study. 
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F.3.4   Wind  Speed  and  Direction  Comparisons 

Stations  were  first  compared  by  examining  the  average  annual 
frequency  of 'occurrence  of  each  wind  direction  (8  point  compass) 
and  calm  conditions.   Each  of  these  nine  values  for  one  station 
in  a  station  pair  was  subtracted  from  the  corresponding  value  for 
the  other  station  to  obtain  their  difference.   The  largest 
difference  across  these  nine  values  was  then  selected  to 
represent  the  maximum  departure  of  each  station  from  the  other  in 
terms  of  wind  direction  frequency.   For  the  two  Toronto  stations, 
the  largest  difference  occurred  for  northerly  winds  which  are 
observed  17  per  cent  of  the  time  at  the  International  Airport  as 
opposed  to  11.7  per  cent  of  the  time  at  the  Downsview  Airport 
for  a  difference  of  5.3  per  cent  (Environment  Canada,  1982).   The 
corresponding  value  for  the  two  Ottawa  stations  was 
3.4  per  cent.   These  values  are  given  in  Exhibit  F.l  along  with 
similar  values  for  all  possible  station  pairs  in  Ontario. 

Wind  speed  characteristics  were  compared  by  examining  average 
annual  wind  speed  values  by  direction  (8  point  compass) .   Again 
the  largest  difference  among  the  eight  wind  speeds  was  selected 
to  represent  the  maximum  departure  of  each  station  from  the  other 
for  each  station  pair.   These  values  are  given  in  Exhibit  F.2  for 
all  Ontario  station  pairs.   For  the  two  Toronto  stations  the 
maximum  wind  speed  difference  was  2.7  km/h,  while  the 
corresponding  value  for  the  Ottawa  stations  was  3.4  km/h. 

F.3.5   Selection  Criteria 

From  these  examples  we  concluded  that  stations  could  be 
considered  to  be  essentially  similar  for  dispersion  modelling 
purposes  if  their  maximum  wind  direction  difference  was  less  than 
5  per  cent  and  if  their  maximum  wind  speed  difference  was  less 
than  3  km/h.   This  rule  is,  if  anything,  too  stringent  since 
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neither  the  Toronto  stations  nor  the  Ottawa  stations  would 
qualify  as  essentially  similar  according  to  our  criteria. 

F.3.6   Selection  of  Locations 

Four  geographical  areas  were  found  to  contain  "essentially 
similar"  stations  that  met  or  almost  met  the  similarity  criteria 
defined  above.   These  areas  included  northwestern  Ontario, 
northeastern  Ontario,  Georgian  Bay  and  southwestern  Ontario. 
Plots  of  average  annual  wind  direction  occurrence  and  average 
annual  wind  speed  by  direction  for  the  candidate  stations  in 
these  regions  are  given  these  regions  in  Exhibit  F.3  through  r.6, 
respectively.   The  two  stations  in  Georgian  Bay  and  Southwestern 
Ontario  obeyed  the  similarity  rule  completely. 

In  northeastern  Ontario,  Timmins  Airport  and  Kapuskasir.g  Airport 
were  essentially  similar  according  to  the  rule  while  the  Timmins 
Airport/Earlton  Airport  pair  (5.3  per  cent  and  4.12  km/h)  and 
the  Earlton  Airport/Kapuskasing  Airport  pair  (3.8  per  cent  and 
3.21  km/h)  did  not.   In  northwestern  Ontario,  Armstrong  Airpcr- 
was  the  station  most  closely  approximating  the  four  others  in 
this  group,  but  only  the  Armstrong  Airport/Sioux  Lookout  pair 
were  essential  similar  according  to  our  selection  rule. 

Despite  the  fact  that  some  northeastern  stations  and  most 
northwestern  stations  had  differences  larger  than  was  specified 
by  the  criteria  of  essential  similarity,  these  stations  are 
located  in  regions  with  very  few  establishments  in  the  economic 
sectors  being  studied.   Hence  the  use  of  additional  dispersion 
patterns  in  those  regions  was  not  justified,  despite  the  slightly 
larger  dissimilarities. 

The  amalgamations  adopted  based  on  the  procedures  described  above 
are  as  follows:   (the  representative  station  is  underlined); 
Toronto  (Toronto  International  Airport.  Toronto  Downsview 


F    -    10 
EXHIBIT   F    -    3 

Wind  Correlations-Northwestern  Ontario 
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Wind   Correlations-Northeastern  Ontario 
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EXHIBIT   F    -    5 
Wind   Correlations-Georgian   Bay 
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EXHIBIT   F    -    6 
Wind   Correlatlons-Southvestern   Ontario 
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Airport) ,  Ottawa  r Ottawa  International  Airport.  Ottawa  Rockclif fe 
Airport) ,  Northwestern  Ontario  r Armstrong  Airport.  Atikokan, 
Geraldton,  PicTcle  Lake,  Sioux  Lookout)  ,  Northeastern  Ontario 
fTimmins  Airrio^t.  Earlton  Airport,  Kapuskasing  Airport)  , 
Georgian  Bay  (Wiarton  Airport,  Gore  Bay  Airport) ,  and 
Southwestern  Ontario  (London  Airport.  Mount  Forest) .   These 
amalgamations  resulted  in  the  removal  of  10  stations  from  the 
data  base. 

Eight  additional  stations  were  also  removed  since  they  are  not 
in  areas  with  significant  sources  of  the  contaminants  studied  in 
this  project.   These  stations  are  Armstrong  Airport,  Big  Trout 
Lake,  Lansdowne  House,  Moosonee,  Muskoka  Airport,  Petawawa 
Airport,  Red  Lake  Airport,  and  Wawa  Airport.   Any  sources 
examined  in  these  regions  use  more  distant  stations  to  define 
average  meteorological  conditions. 

The  stations  retained,  which  together  are  taken  to  define  the 
meteorological  conditions  for  the  province,  are  Kenora  Airport, 
Kingston  Airport,  London  Airport,  North  Bay  Airport,  Ottawa 
International  Airport,  Sault  Ste.  Marie  Airport,  Simcoe,  Sudbury 
Airport,  Thunder  Bay  Airport,  Timmins  Airport,  Toronto 
International  Airport,  Trenton  Airport,  Wiarton  Airport  and 
Windsor  Airport. 


F.4   COMPILATION  OF  METEOROLOGICAL  DATA 


F.4.1  Need  to  Select  a  Representative  Year 

While  the  rationalization  of  observing  stations  to  be  included  in 
the  analysis  reduced  the  data  to  be  analyzed  by  more  than  50  per 
cent,  the  data  burden  was  still  too  large  for  the  project.   To 
accomplish  a  further  reduction  it  was  decided  to  replace  the 
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30  year  period  of  record  used  to  characterize  meteorological 
conditions  with  a  representative  year  that  closely  approximates 
the  30  year  petiod.   This  representative  year  yields  an  annual 
dispersion  field  very  close  to  that  obtained  from  the  3  0  year 
history  for  each  of  the  14  stations. 

For  the  purposes  of  the  project  only  annual  statistics 
summarizing  the  dispersion  pattern  are  required,  the  large  number 
of  hourly  observations  in  one  year  ensure  that  statistics  derived 
from  these  hourly  observations  do  not  vary  greatly  from  year  to 
year.   Consequently,  any  year  in  the  data  base  could  be  used  to 
represent  typical  conditions  if  all  input  data  were  hourly 
observations.   Such,  however,  is  not  the  case. 

F.4.2   Snow  Cover 

Specifically  the  presence  or  absence  of  snow  on  the  ground  is 
not  subject  to  the  law  of  large  numbers  in  the  way  that  hourly 
observations  are.   The  period  of  time  that  snow  is  present  for 
any  station  varies  a  great  deal  from  year  to  year.   In  the 
Ministry's  model  the  presence  of  snow  on  the  ground  is  assumed 
to  reflect  all  short  wave  radiation  back  into  space,  with  the 
result  that  convective  conditions  cannot  occur  (Ontario  Ministry 
of  the  Environment,  1987) .   The  disallowance  of  convective 
conditions  for  an  unusually  long  or  unusually  short  period 
relative  to  the  30  year  average  may  well  result  in  a 
significantly  different  dispersion  field.   Accordingly,  it  was 
decided  to  seek  a  single  year  which,  for  all  14  stations,  the 
winter  snow  cover  period  was  as  close  as  possible  to  the  average 
obtained  using  the  entire  period  of  record  for  the  station. 

The  time  period  for  which  snow  was  on  the  ground  for  at  least 
seven  consecutive  days  was  determined  for  each  of  the  14 
stations  for  all  years  of  each  stations 's  climatological  record. 
Significant  year  to  year  variability  was  found.   For  example,  in 


F  -  16 

Kenora  a  30  year  data  base  indicated  that  on  average  seven 
consecutive  days  of  snow  cover  began  on  November  11  and  ended  on 
April  11  for  a ' total  of  152  days.   However,  in  1968  snow  was 
present  for  only  126-  days  while  in  1965  snow  was  present  for  179 
days.   Other  stations  showed  similar  variability  in  terms  of  the 
length  of  the  snow  cover  season. 

F.4.3   Selection  of  the  Representative  Year 

Only  one  year,  198  6,  had  the  time  with  snow  on  the  ground  within 

one  standard  deviation  of  the  long  term  average  for  all 

14  stations.   Thus  198  6  was  selected  as  the  representative  year 

for  the  calculation  of  average  annual  dispersion  statistics. 

F.4.4   Compilation  of  Climatoloqical  Input  Data;   Hourly  Station 
Data 

Hourly  meteorological  observations  for  the  representative  year 
for  each  of  the  14  centres  (totalling  500,000  observations)  were 
purchased  from  the  Canadian  Climate  Centre,  Atmospheric 
Environment  Service  and  formatted  for  input  to  the  Ministry's 
METCO  program. 

The  METCO  program,  provided  by  the  Ministry,  allows  the  use  of 
data  formatted  for  use  in  the  single  source  GAS  model  in  the 
newer  multi-source,  multi-receptor  version  of  the  model.   The 
program  also  allows  disabling  of  upper  air  wind  speeds  where 
desirable. 

This  program  was  used  as  the  basis  for  a  general  conversion 
program  (METCOR)  written  by  the  study  group  to  automate  the 
interface  of  Environment  Canada  Surface  Data,  U.S.  and  Canadian 
Upper  Air  Data,  the  snow  record,  data  decision-making,  flagging 
and  reformatting  of  input  data.   The  METCOR  program  produced 
14  complete  records  for  the  representative  year,  ready  for  input 
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to  the  Ministry's  MET  program  (meteorological  parameter 
calculation  program) .   These  14  records  represent  average 
meteorological' conditions  for  the  Province  of  Ontario  for  the 
period  1957  td '1986  inclusive.   The  MET  program  generates  the 
meteorological  parameters  required  as  input  by  the  dispersion 
(GAS)  model. 

F.4.5   GMT  Conversion 

Canadian  surface  weather  data  are  recorded  in  local  standard  time 
for  each  station.   A  program  segment  was  written  to  convert  local 
standard  time  (date)  to  Greenwich  Mean  Time  (GMT)  for  input  to 
the  GAS  model. 

F.4.6  Cloud  Cover 

Canadian  surface  weather  data  expresses  cloud  cover  and  cloud 
opacity  in  lOths.   Cloud  cover  and  cloud  opacity  affect  the 
occurrence  of  convective  conditions.   Cloud  opacity,  where 
available,  is  the  better  determinator  of  convective  conditions 
and  was  used  in  this  study.   A  program  segment  was  used  to 
convert  from  lOths  to  8ths  for  input  to  the  GAS  model. 

F.4.7   The  Snow  Record 

The  most  significant  determinators  of  convective  conditions  are 
the  presence  or  absence  of  snow  cover  and  the  presence  or 
absence  of  daylight.   Most  winters  are  characterized  by 
alternating  periods  of  "snow-on"  and  "snow-off".   Wide  swings 
occur,  throughout  the  30-year  record,  of  starting  and  ending 
dates  -  dates  of  snow  cover,  and  dates  of  bare  ground.   The  snow 
cover  period  may  also  be  interrupted  by  periods  of  bare  ground. 

The  30-year  record  for  each  centre  was  reduced  to  a  mean  start 
and  end  date  of  snow  cover.   The  definition  of  snow  cover  used  is 
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the  period  from  which  persistent  (7-day  duration) ,  2  cm  snow 
cover  was  first  recorded  to  its  last  recorded  date.   This 
definition  had  the  effect  of  reducing  the  length  of  the  snow 
season.   The' reduction  in  the  length  of  the  season  offsets,  to 
some  degree,  the  increase  in  the  number  of  "snow-on"  days  which 
result  from  calling  the  snow  cover  period  "continuous". 

The  3  0-year  mean  period  was  used  as  the  "snow-on"  period  for  the 
representative  year.   A  program  segment  was  written  to  automate 
snow  cover  distribution  in  the  model. 

F.4.8  Wind 

Canadian  surface  wind  speed  is  in  units  of  km/hr.   A  program 
segment  was  written  to  convert  this  to  m/sec  for  the  GAS  model. 
Calm  conditions  were  assigned  a  wind  speed  value  of  0  m/sec  where 
no  value  was  provided.   The  MET  and  GAS  modules  increase  this  to 
a  1.5  m/sec  minimum  in  appropriate  cases.   Upper  air  wind  speeds 
were  disabled  in  the  model. 

F . 4 . 9   Wind  Direction 

The  representative  year  chosen,  being  recent,  was  available  in 
10  degree  compass  resolution  and  was  thus  compatible  with  the 
GAS  model  interpolation  routine  (wind  direction  +5  degree 
randomization) .   A  program  segment  was  written  to  multiply  the 
AES  compass  (0  to  36)  by  10.   In  Environment  Canada  format,  36  is 
North  and  0  is  calm.   A  program  segment  was  written  to  reassign 
36  to  0  for  the  GAS  model.   Initial  0  (calm)  values  were 
assigned  the  first  previous  non-calm  wind  direction. 

F.4.10  Upper  Air  Data 

No  Canadian  upper  air  stations  exist  in  proximity  to  most 
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Ontario  emission  centres.   Thus,  it  was  necessary  to  obtain  U.S. 
near-border  data  to  use  with  Ontario  surface  data. 

F.4.11  Mixinct^ Heights;   U.S.  Locations 

Mixing  heights  were  purchased  for  the  representative  year  from 
the  U.S.  National  Climatic  Data  Centre  in  North  Carolina.   The 
U.S.  computer  program  used,  generates  a.m.  and  p.m.  mixing 
heights  using  upper  air  soundings  combined  with  surface  data  for 
the  same  station  (Doty,  1983).   The  normal  U.S.  practice  is  to 
combine  surface  data  from  the  surface  location  of  interest. 
However,  since  the  surface  locations  of  interest  are  all 
Canadian  centres,  this  was  not  possible.   The  use  of  matching 
upper  air/surface  data  for  the  same  station  was  considered  a 
reasonable  compromise. 

Twice  daily,  mixing  heights  for  the  entire  year  were  obtained 
for  Flint,  Michigan;   Sault  Ste.  Marie,  Michigan;   International 
Falls,  Minnesota;   and  Buffalo,  New  York. 

F.4,12  Mixing  Heights;   Canadian  Locations 

Since  mixing  heights  cannot  be  purchased  from  Environment  Canada 
for  Canadian  Upper  Air  Stations,  raw  soundings  data  were 
purchased.   The  only  Canadian  centre  in  proximity  to  the 
locations  to  be  modelled  is  Maniwaki,  Quebec. 

Since  Maniwaki,  Quebec  surface  data  were  not  purchased,  Ottawa 
surface  data  were  used  with  the  Maniwaki  upper  air  data  to 
generate  mixing  heights.   This  is  in  accordance  with  U.S. 
practice.   A  program  was  written  to  adjust  base  temperatures  for 
base  elevations,  determine  maximum  surface  temperatures  for  each 
period,  and  calculate  the  elevation  of  intersection  between  the 
sounding  temperature  profile  and  the  dry  adiabatic  lapse  curve 
for  each  surface  temperature  maximum. 
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A  linear  interpolation  routine  was  written  to  refine 
intersections  Which  occurred  between  sounding  elevations. 
Negative  eleVations  were  assigned  a  value  of  zero,  which  was 
flagged  by  the  revised  MET  program  for  cases  which  would 
ordinarily  be  convective.   In  such  cases,  a  zero  mixing  height 
resulted  in  an  invalid  MET  case  in  the  GAS  model  execution 
(i.e.,  no  calculation). 

Ottawa  was  chosen  as  the  surface  centre  because  it  is 
geographically  central  to  the  area  of  Ontario  for  which 
Maniwaki,  Quebec  upper  air  data  would  be  deemed  representative 
(i.e.,  Timmins,  North  Bay,  Ottawa,  Kingston).   Time  and  budget 
did  not  permit  a  centre-by-centre  match  of  local  surface  and 
Maniwaki  upper  air  data. 

F.4.13  Missing  Data  Flags 

In  the  case  of  missing  meteorological  data  a  flag  value  of  IDATA 
=  999  was  assigned.   This  generates  an  invalid  MET  case  for  the 
GAS  run  (i.e.,  no  calculation  for  that  case). 

F.4.14  Negative  or  Missing  Mixing  Height  Flags 

In  cases  of  negative  or  missing  mixing  heights,  IDATA  was 
assigned  a  value  of  987.   A  short  addition  to  the  MET  program  was 
written  such  that  if  IDATA  =  987  and  meteorological  conditions 
were  convective,  IFLAG  would  be  assigned  a  value  of  999  (an 
invalid  case) . 

F.4.15  Limitations  of  the  Meteorological  Database 

Some  centres,  despite  catalogue  indications,  were  found  to  have 
incomplete  meteorological  records.   For  example,  in  Simcoe, 
measurements  were  only  available  for  weekday  periods.   Results  of 
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modelling  with  such  data  are  less  reliable  than  results  for  more 
detailed  input. 

The  use  of  Ufi^'er  Air  Data  from  geographically  removed  centres  is 
a  somewhat  dubious  procedure.   However  it  is  one  that  is  without 
alternative  since  upper  air  stations  are  so  rare.   For  example, 
using  Buffalo  upper  air  data  for  Toronto  leads  to  questions  of 
how  representative  upper  air  data  for  a  centre  on  the  east  end  of 
a  lake  is  of  another  centre  on  the  north  shore  of  a  different 
lake.   Errors  would  occur  only  in  the  case  of  convective 
conditions  and  are  estimated  to  yield  an  output  variance,  on  an 
annual  average,  of  no  more  than  5  to  10  per  cent.   The  best 
centre,  by  this  criterion,  is  Sault  Ste.  Marie,  Ontario,  since 
the  Sault  Ste.  Marie,  Michigan  upper  air  station  is  so  close  to 
the  Canadian  surface  weather  station. 

There  are  many  surface  stations  in  northwestern  Ontario.   Kenora 
was  chosen  to  represent  this  area  based  on  its  population 
density.   It  is  estimated  that  this  choice  would  affect  the 
dispersion  of  pollutants  in  other  centres  (e.g.,  Armstrong)  to  a 
smaller  degree  (5  to  10  per  cent) .   It  was  not  expected  that 
large  concentration  excursions  would  occur  and  none  were 
observed. 

The  dispersion  model  (GAS)  does  not  deal  with  calm  conditions  "as 
are"  but  rather  increases  wind  speed  to  1.5  m/sec.  for  applicable 
cases.   In  addition,  Environment  Canada  merely  flags  wind 
direction  in  calm  cases  (since  by  definition  of  calm,  there  is  no 
wind  to  have  a  direction) .   To  mesh  what  the  model  does,  with 
what  the  source  data  doesn't  provide,  a  wind  direction  was 
assigned  in  each  case  of  calm.   The  choice  of  direction  was  that 
of  the  first  previous  "non-calm"  wind  direction.   This  perturbs 
the  short  term  wind  direction  statistic  less  than  introduction  of 
an  artificial  randomization  function.   The  frequency  of 
occurrence  of  calm  conditions  is  estimated  at  less  than 
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2  per  cent,  so  the  impact  of  this  assumption  on  results 
considered  to  be  small. 

..  .  -''*'■''  • 

F.5   EMISSION  SOURCE  CHARACTERISTICS 


The  GAS  model  requires,  in  addition  to  climatological  data,  input 
with  respect  to  emission  source  characteristics.   The 
characteristics  of  the  several  thousand  emissions  sources  are, 
with  very  rare  exceptions,  not  known.   Hence  these 
characteristics  must  be  estimated.   Since  each  of  the 
14  dispersion  patterns  could  be  applied  to  several  hundred 
sources,  the  estimated  characteristics  should  be  typical  of  the 
variety  of  industrial  establishments  emitting  the  contaminants. 

F.5.1  Building  Sources 

Of  the  several  thousand  commercial  emissions  sources  in  Ontario 
subject  to  the  proposed  revisions  to  Regulation  3  08,  by  far  the 
most  numerous  are  building  sources.   Such  sources  are  generally 
defined  as  sources  from  which  fugitive  emissions  emanate  or  stack 
emissions  emerge  at  elevations  less  than  twice  the  height  of  the 
building. 

To  precisely  model  emissions  from  such  sources  the  exact 
dimensions  of  the  building  and  the  stack  must  be  known.   In  cases 
of  complex  building  shape  (i.e.,  non-rectangular)  or  complex 
vicinity  (i.e.,  larger  buildings  or  variable  terrain) 
mathematical  modelling  is  inappropriate  and  physical  modelling  is 
preferred.   Physical  modelling  is  not  possible  for  this  analysis. 
Hence  it  is  necessary  to  assume  that  all  buildings  are  isolated, 
rectangular  buildings  with  a  stack  less  than  twice  the  height  of 
the  building. 
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For  this  study,  no  information  was  available  on  building  or 
building-stack  dimensions  for  any  source.   To  estimate  ambient 
concentration,'  each  source  was  assumed  to  be  a  building  of  7  m 
height  and  100', m  width  with  a  stack  of  less  than  14  m  elevation, 
emitting  pollutants  at  0  velocity  and  at  ambient  exit 
temperature.   Such  a  building  is  a  typical  "industrial  park" 
source.   This  assumption  leads  to  very  high  ambient 
concentration  uncertainties  (i.e.,  order  of  magnitude)  within 
500  m  of  the  source  for  source  buildings  of  other  configuration 
but  is  a  far  better  source  assumption  than  a  zero  elevation 
emitter.   Reasonable  concentration  estimates  beyond  500  m  are 
achieved. 

Since  no  data  were  available  on  individual  elevated  and  buoyant 
(exit  temperature  above  the  ambient  temperature)  elevated  source 
characteristics,  such  sources  have  been  represented  as  building 
sources.   The  limitations  of  this  approach  are  discussed  in  the 
next  section. 

F.5.2   Implications  of  the  Assumed  Building  Source 
Characteristics 

If  an  emission  is  ground-based,  the  building  source  assumption 
will  yield  artificially  high  initial  dilution  within  500  m  of  the 
origin.   If  an  emission  is  elevated,  and  possibly  buoyant,  the 
building  source  assumption  will  yield  artificially  high 
concentrations  within  500  m  of  the  source.   If  the  source  is  less 
than  about  30  m  in  elevation  and  of  moderate  buoyancy,  beyond 
500  m  concentrations  will  generally  fall  within  the  same  order  of 
magnitude  with  the  highest  discrepancies  found  at  the  perimeter 
(near  the  zero  concentration  area).   In  general,  the  higher  and 
hotter  the  source,  the  larger  the  discrepancies. 

Distribution  inaccuracies  resulting  from  the  assumptions  can  be 
illustrated  by  the  fact  that  the  worst  case  meteorology  for  an 
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elevated  source  occurs  under  convective  atmospheric  conditions, 
whereas  the  worst  case  meteorology  for  a  building  source  occurs 
under  stable  atmospheric  conditions.   This  type  of  possible 
inaccuracy  contributed  to  the  decision  to  use  concentration 
averages  of  8,760  hourly  meteorological  cases  rather  than 
averages  over  shorter  time  frames.   Where  sources  approached  the 
characteristics  of  the  "typical  building  source"  the  accuracy  of 
the  annual  average  would  be  unassailable.   In  other  cases,  the 
annual  average  would  have  a  "smoothing"  effect  where  periods  of 
larger  differences  would  be  offset  to  some  degree  by  periods  of 
smaller  differences  yielding  a  better  estimate  of  average 
concentration.   This  minimizes  the  limitations  of  these  source 
assumptions. 

F.5.3   Elevated  Sources 

Elevated  sources,  by  definition,  emit  beyond  twice  the  elevation 

of  the  highest  part  of  an  attached  or  adjacent  building. 
Depending  on  the  stack  exit  temperature  the  emission  plume  may  be 
more  or  less  buoyant  than  the  surrounding  air.   Accurate 
modelling  requires  that  all  of  the  stack  parameters  be  known. 
These  parameters  are:   stack  height,  exit  diameter,  exit 
velocity,  exit  temperature  (and  emission  rate) . 

Such  information  was  not  collected  for  any  source  in  Ontario  as 
part  of  this  study.   In  addition,  the  use  the  proposed  new 
Ontario  models  entirely  precluded  use  of  individual  source 
information  of  any  kind.   If  one  took  2,000  sources  for  example, 
and  ran  one  years'  meteorology,  using  an  IBM-AT  with 
math-coprocessor,  3  6,000  hours  of  computer  time  would  have  been 
required  (i.e.,  4.1  years). 
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F.5.4   Near-Lake  Sources 

Discussions  wit;h  Ministry  staff  led  to  consideration  that 
lakeshore  effects  should  be  dealt  with  in  modelling  certain 
sources. 

Predominant  sources  to  which  lakeshore  effects  would  most  often 
be  encountered  were  identified  as  Lakeview  Generating  Station  in 
Mississauga,  and  Nanticoke  Generating  Station.   The  present 
development  version  of  the  Ministry's  proposed  Lakeshore  Model 
has  not  been  integrated  with  the  general  GAS  model.   Such 
integration  is  necessary  because  a  lakeshore  effect  can  only 
occur  under  specific  conditions  which  will  only  arise  at  certain 
times  during  the  year.   At  other  times  dispersion  would  occur  as 
per  the  GAS  model. 

We  endeavoured  to  perform  the  model  integration;   however,  the 
project  budget  was  exceeded  as  a  result  and  this  work  remains 
incomplete.   Hence,  lakeshore  effects  are  not  included  in  the 
output  of  this  study. 

The  error  involved  in  ignoring  lakeshore  effects  for  these  two 
sources  would  be  an  under-estimation  in  the  range  of  1  to  3  times 
of  annual  average  concentrations  within  3  km  onshore  from  these 
two  sources. 


F.6   THE  RECEPTOR  GRID 


The  Ministry's  grid-generator  program  RECGEN  was  used  to  generate 
a  receptor  grid  network  from  0  to  24  km  of  2  km  resolution.   A 
finer  grid  of  500  m  and  1,000  m  grid  points  was  superimposed  on 
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the  coarse  grid  to  permit  better  resolution  in  the  near-source 
area  where  most  of  the  impact  of  building  sources  is  felt. 

The  GAS  model '■fcalculated  running  average  (annual  average  in  this 
study)  concentrations  for  each  receptor  point  (grid  coordinate) . 
A  surface  roughness  of  1.0  m  and  a  receptor  elevation  of  0.0  m 
was  used  in  all  cases. 

Due  to  the  assumptions  relating  to  source  characteristics  the 
resulting  ambient  concentrations  within  500  m  are  most 
uncertain.   As  an  approximation,  where  receptor  grids 
intersect/ overlap  (for  different  sources), if  the  source  is  known 
to  be  elevated  the  within-500  m  area  could  be  taken  to  be  of  zero 
concentration.   Where  the  source  was  not  known  to  be  elevated, 
the  within-500  m  area  could  be  taken  to  be  of  concentration  equal 
to  that  at  the  500  m  boundary.   The  alternative  would  be  to 
ignore  all  areas  of  overlap  within  500  m  of  the  source  (which  in 
the  case  of  Toronto  might  mean  ignoring  most  of  the  city  for 
some  pollutants) . 

Fourteen  dispersion  patterns  were  generated;   one  for  each 
location  selected.   The  source  characteristics  are  assumed  to  be 
the  same  for  all  14  patterns.   An  emission  rate  of  1  gram/sec  is 
used  in  each  case.   Since  the  model  is  linear  with  respect  to 
dispersion,  the  calculated  concentration  at  each  grid  coordinate 
is  scaled  by  the  emission  rate  for  the  specific  source.   For 
example,  if  the  concentration  at  twelve  kilometers  due  east  of 
the  source  is  calculated  by  the  model  to  be  0.0158  micrograms/ 
cu.m.  and  a  given  source  emits  a  contaminant  at  a  rate  of 
5  grams/sec  then  the  ambient  concentration  twelve  km  due  east  of 
that  source  will  be  0.0790  micrograms/cu.m. 
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F.7   QUALITY  ASSURANCE 


F.7.1  IriTsut- 'Quality  Control 

Considerable  processing,  conversion,  and  flagging  of  data  were 
required  prior  to  running  the  dispersion  model.   Conversion, 
processing,  and  flagging  programs  were  all  validated  manually 
using  Toronto  and  Simcoe  meteorological  data.   MET  program 
output  (GAS  program  input)  was  reviewed  in  detail  (i.e.,  for  the 
entire  one  year  hourly  record)  prior  to  GAS  program  execution. 
Byte  volumes  of  records  for  other  centres  were  compared  for 
consistency  and  random  hourly  records  were  inspected  for  other 
centres. 

All  centres  used  the  same  input  files,  except  for  meteorological 
data  where  individual  centre  files  were  used.   To  maintain  file 
control,  all  non-.MET  files  (i.e.,  .RUN,  .FUN,  .REC,  .SRC  and 
blank  .AVG  files)  were  located  in  a  write-protected  directory.   A 
DOS  batch-program  was  written  which  transferred  these  files 
(renamed)  to  the  execution  directory,  prompted  for  the  correct 
•MET  file,  verified  the  presence  of  these  files,  and  executed  the 
GAS  model  program. 

F.7. 2   Output  Quality  Control 

The  number  of  valid  and  invalid  (uncalculated)  MET  cases  was 
inspected  on  completion  of  each  run  and  compared  with  the  prior 
knowledge  of  the  relative  occurrence  of  missing  and  flagging 
data.   Exhibit  F.7,  the  Run  Record,  illustrates  this. 

All  output  (.AVG)  files  were  imported  to  a  LOTUS  Symphony  file, 
compressed,  and  compared.   The  coordinates  of  the  maximum 
concentration  receptor  were  inspected  in  each  case. 
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The  relative  concentrations  of  the  maxima  were  compared.  With 
each  centre  using  the  same  building  source  and  unit  emission 
rate,  the  maxilna  should  be  within  the  same  order  of  magnitude. 
In  addition  thl©  concentrations  at  every  grid  point  were  summed 
for  each  centre  and  the  sums  compared  between  centres .  This  test 
on  distribution  should  also  yield  results  in  the  same  order  of 
magnitude . 

All  of  the  above  tests  were  passed.   Exhibit  F.8  presents  the 
first  page  of  this  summary  file  as  an  example. 

A  DOS  batch  program  was  written  to  copy  the  appropriate  files  to 
diskette  for  back-up  and  to  save  output  files  to  the  appropriate 
directory.   The  DOS  programs  ensured  that  file  integrity  was 
maintained  during  movement  and  handling. 
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G . 1   INTRODUCTION 


A  vital  tool  in  the  analysis  of  the  public  benefits  of  proposed 
revisions  to- Regulation  3  08  used  was  DPA's  spatial  analysis 
system  (SPANS) .   SPANS  is  a  geographic  information  system. 
Different  geographic  information  systems  specialise  in  performing 
different  tasks;   some  work  best  at  creating  cartographic  maps. 
SPANS'  specialty  is  analysis  and  manipulation  of  spatial  data; 
data  which  relate  to  points  or  areas  on  maps.   It  was  ideally 
suited  to  this  assignment  which  involved  locating  sources  of 
emissions,  the  dispersion  of  those  emissions,  and  identification 
of  the  population  and  environmental  resources  affected. 

This  appendix  describes  some  of  the  features  of  SPANS  and  the  way 
it  was  used  in  the  project. 


G.2   WHAT  IS  SPATIAL  ANALYSIS? 


Spatial  analysis  is  the  ability  to  investigate  the  spatial 
dimension  of  a  problem.   Many  problems  have  a  spatial  dimension. 
For  instance,  the  analysis  of  the  public  benefits  of  reduced  air 
emissions  depend  upon  where  contaminants  are  deposited  and  the 
concentration  at  each  point  because  that  determines  the  resources 
affected.   Traditionally,  the  spatial  dimensions  of  problems  are 
addressed  by  making  assumptions  (e.g.,  trade  area  boundaries)  or 
using  manual  procedures  (map  overlays)  that  limit  the  number  of 
options  considered. 

DPA's  computerized  spatial  analysis  system  allows  the  spatial 
components  of  a  problem  to  be  addressed  in  much  more 
sophisticated  and  convenient  ways.   To  understand  how  the 
spatial  analysis  system  works  some  background  information  on 
geographic  infoirmation  systems  is  useful. 
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G.3   GEOGRAPHIC  INFORMATION  SYSTEMS 


Computerized  geographic  information  systems  (GIS)  are  designed  to 
input,  store,  retrieve,  manipulate  and  display  spatial 
information.   They  differ  primarily  in  the  manner  in  which  they 
store  data  and  the  kind  of  manipulations  they  perform  best.   In 
turn  these  capabilities  determine  the  size  of  computer  needed. 

G.3.1   Data  Storage 

A  spatial  feature  has  a  location,  size  and  shape.   The  size  and 
shape  can  range  from  a  point  (e.g.,  a  meterological  station)  to 
an  irregular  area  (e.g.,  a  watershed).   In  addition  to  these 
spatial  characteristics,  the  GIS  must  store  a  range  of  attributes 
associated  with  the  site:   average  annual  temperature,  soil  type, 
slope,  population,  suitability  for  swimming,  etc.   Clearly,  the 
data  storage  demands  can  quickly  become  very  large,  creating  a 
need  for  larger  computers. 

The  most  common  approach  to  accurately  specify  spatial  features 
is  to  define  their  boundaries.   A  series  of  line  segments  are 
linked  to  form  the  boundary  of  an  area.   Obviously,  the  more 
irregular  the  area  the  larger  the  storage  requirements.   Some 
analytical  processes,  such  as  map  overlays,  are  very  complex 
with  these  systems.   Line  segments  must  be  used  to  generate  the 
areas,  the  areas  are  overlaid,  the  portions  that  overlap  are 
determined  and  a  new  set  of  boundary  line  segments  is  generated 
for  storage.   These  systems  typically  require  mini-computers. 

The  usual  approach  to  minimising  storage  needs  is  to  pre-define  a 
fixed  set  of  areas.   Typically  this  is  a  grid  (raster)  system, 
although  other  systems,  such  as  census  subdivisions,  can  be  used. 
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The  fixed  areas  and  locations  defined  by  such  systems  yield  major 
benefits  in  terms  of  data  storage.   The  cost  is  loss  of  accuracy. 
Watersheds,  slopes,  soils,  rainfall  and  many  other  attributes  do 
not  vary  acctJrtiing  to  a  fixed  grid  pattern. 

SPANS  uses  a  quadtree  structure  to  store  data.   This  offers  a 
combination  of  the  best  features  of  a  boundary-based  system  and 
grid  system.   The  quadtree  structure  divides  the  total  area  into 
progressively  smaller  quadrants,  until  a  quadrant  has  uniform 
characteristics.   This  results  in  area  units  of  variable  size,  so 
a  large  area  with  uniform  attributes  is  stored  as  such,  while 
much  smaller  area  units  are  used  to  accurately  define  the 
boundaries  of  spatial  features,  such  as  forest  cover,  lakes,  and 
built  up  area.   The  overall  data  storage  requirements  are 
similar  to  those  for  a  fixed  area  system,  but  the  level  of 
accuracy  is  much  greater.   The  system  operates  on  a  micro 
computer. 

G.3.2   Manipulation  Capabilities 

The  manipulation  capabilities  of  a  GIS  system  reflect  the  needs 
of  users  and,  to  some  extent,  the  data  storage  techniques 
employed.   The  available  systems  can  be  broadly  defined  as 
cartographic  or  thematic  in  orientation.   The  principal  uses  of 
cartographic  systems  are  to  record,  process  revisions  to,  and 
prepare  updated  maps.   Although  these  systems  have  some 
analytical  capabilities,  such  as  the  ability  to  perform  overlays, 
the  emphasis  is  on  cartographic  properties  such  as  accuracy, 
labelling,  changes  to  scale  and  projection,  etc. 

Conversely  thematic  systems  emphasize  analytical  capabilities  but 

at  the  cost  of  some  cartographic  ability  such  as  labelling. 

SPANS  was  specifically  designed  as  an  analytical  system.   Its 
capabilities  include: 
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.  Easy  access  to  a  variety  of  input  data; 

.  Extensive  data  manipulation  capabilities; 

.  Complex  overlays  of  multiple  maps; 

.  Proximity  mapping; 

.  Area  analysis  reports; 

.  Use  of  satellite  imagery;   and 

.  Diverse  data  display  capabilities. 

SPANS  grew  out  of  needs  encountered  in  the  course  of  DPA ' s 
consulting  work  over  the  past  decade.   The  system  was  developed, 
refined  and  enhanced  over  many  years  and  a  wide  variety  of 
applications. 


G.4   USE  OF  SPANS  IN  THIS  PROJECT 


A  critical  requirement  in  the  calculation  of  the  public  benefits 
of  reduced  air  emissions  is  the  ability  to  estimate  the  change 
in  concentration  of  a  given  contaminant  at  any  point  and  to 
determine  the  environmental  resources  that  will  benefit  from  this 
reduced  exposure.   SPANS  provides  this  capability. 

G.4.1   Location  of  Emission  Sources 

SPANS  is  used  to  locate  each  emission  source.   The  list  of 
establishments  was  compiled  from  the  Made  in  Ontario  data  base 
maintained  by  the  Ministry  of  Industry,  Trade  and  Technology  and 
Statistics  Canada's  list  of  establishments.   The  Made  in  Ontario 
data  base  includes  the  postal  code  and  these  establishments  were 
located  at  the  (six  digit)  postal  code  centroid.   The  centroid 
locations  for  postal  codes  were  purchased  from  Statistics 
Canada.   The  Statistics  Canada  list  only  gives  the  municipality. 
These  establishments  were  located  at  the  centroid  for  the 
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municipal  unit.   The  municipality  centroid  was  calculated  using 
SPANS. 

*"  "i.  1  .  ,  , 

To  facilitate  the  calculation  of  depositions  from  different 

sources,  a  standard  4  km  grid  was  created.   Each  establishment 

was  moved  to  the  nearest  grid  point,  a  distance  of  less  than 

2  km.   Given  that  the  initial  location  was  the  postal  code  or 

municipality  centroid,  this  did  not  introduce  a  significant  bias 

or  error  into  the  analysis.   Grid  points  with  numerous  (over  10) 

establishments  were  examined.   The  larger  firms  assigned  to  that 

grid  point  were  contacted  and,  if  appropriate,  shifted  to  other 

grid  points.   However,  it  remains  true  that  establishments  are 

artificially  concentrated  on  some  grid  points  due  to  the  lack  of 

accuracy  of  the  address  information.   The  establishments  from  the 

Statistics  Canada  list  are  concentrated  at  fewer  locations 

because  there  are  fewer  municipalities  than  postal  codes. 

The  attributes  stored  for  each  establishment  are  its: 
.  Canadian  SIC  code; 
.  estimated  employment  (or  other  size  characteristic) ; 

.  estimated  current  emissions  by  contaminant; 

.  estimated  emissions  by  contaminant  for  each  implementation 

scenario;   and 
.  an  indicator  to  identify  the  appropriate  dispersion  pattern. 

G.4.2   Determination  of  Concentrations 

The  dispersion  modelling  (see  Appendix  F)  calculated  14 
dispersion  patterns  for  different  parts  of  the  province.   The 
dispersion  model  assumed  an  emission  rate  of  1  gm/sec  and 
calculated  the  resulting  concentration  (micrograms/cu. m)  on  a 
2  km  grid  covering  an  area  24  km  in  all  directions  from  the 
source. 
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To  keep  the  calculations  manageable  —  65  contaminants, 
over  3,500  establishments,  5  scenarios  —  it  was  necessary  to 
adopt  a  coarser  grid.   Thus,  a  4  km  grid  rather  than  a  2  km  grid 
was  used.   Sln'ce  the  sources  had  been  moved  to  points  on  a  4  km 
grid,  the  concentration  estimates  also  corresponded  to  points  on 
this  grid. 

The  analysis  proceeds  on  a  contaminant  by  contaminant  basis.   If 
an  establishment  is  a  source  of  the  contaminant,  the  emissions 
for  the  specified  scenario  are  converted  to  an  emission  rate. 
The  appropriate  dispersion  pattern  is  selected.   The  origin  of 
the  dispersion  pattern  is  adjusted  to  the  source  location.   Every 
point  in  the  dispersion  grid  is  similarly  adjusted  to  points  on 
the  province-wide  grid.   The  concentration  at  each  point  on  the 
dispersion  grid  is  scaled  by  the  emission  rate  and  added  to  the 
existing  concentration  at  that  point  on  the  province-wide  grid. 

When  this  process  has  been  completed  the  grid  points  reflect  the 
estimated  concentration  of  the  contaminant  at  that  point  due  to 
emissions  from  all  establishments  under  the  given  scenario.   The 
concentration  at  each  grid  point  is  assumed  to  apply  to  a  4  km 
square  area  with  the  grid  point  at  its  centre. 

G.4.3   Calculation  of  Ambient  Concentration 

The  dispersion  model  estimates  the  concentration  of  a 
contaminant  due  to  a  specific  source.   Many  of  the  contaminants 
analysed  are  also  emitted  by  sources  other  than  those  considered 
in  this  study.   Thus,  the  ambient  concentration  due  to  all 
sources  is  likely  to  be  higher  than  the  estimated  emissions  from 
the  Regulation  3  08  sources  covered  by  this  study. 

Data  on  total  ambient  concentrations  are  collected  by  monitoring 
stations.  Two  networks  of  monitoring  stations  exist  in  Ontario. 
One  network  has  stations  located  in  areas  where  air  quality 
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problems  are  believed  to  exist.   The  other  network  has  stations 
located  more  uniformly  across  the  province  to  collect  data  to 
monitor  acidic!  deposition.   Each  monitoring  station  collects 
information  tin',  ambient  concentrations  of  a  number  of 
contaminants.   Of  the  96  contaminants,  monitoring  data  are 
available  for  cadmium,  chromium,  lead,  manganese,  nickel, 
nitrogen  oxides,  nitrous  oxide,  sulphur  dioxide,  total 
particulate  matter  (PM)  and  total  VOC's. 

For  each  of  these  contaminants  except  total  VOC's,  the  ambient 
concentration  as  recorded  by  the  monitoring  station  was  compared 
with  the  concentration  at  the  monitoring  station  location 
estimated  using  the  dispersion  model.   The  result  is  a  series  of 
"residuals"  (one  for  each  monitoring  station)  that  represent  the 
ambient  concentration  due  to  non-Regulation  3  08  sources.   These 
residuals  were  examined  individually  for  each  contaminant.   Based 
on  this  examination,  one  (or  more)  values  of  the  ambient 
concentration  due  to  non-Regulation  308  sources  were  established. 
The  value (s)  and  the  rationale  for  the  selection  is  presented  in 
Chapter  4  in  the  section  relating  to  the  specific  monitored 
contaminant. 

The  value (s)  of  the  concentration  due  to  non-Regulation  308 
sources  is  added  to  the  estimated  concentration  due  to 
Regulation  308  sources  at  each  grid  point  for  all  scenarios. 
This  converts  the  value  at  each  grid  point  from  an  estimate  of 
the  concentration  due  to  the  specific  Regulation  308  sources 
analysed  to  an  estimate  of  total  ambient  concentration.  And  it 
permits  an  analysis  of  the  impact  of  the  proposed  changes  to 
Regulation  308  on  total  ambient  concentration  of  the  given 
contaminant.  The  ambient  concentration  at  each  grid  point  is 
assumed  to  apply  to  a  4  km  square  area  with  the  grid  point  at 
its  centre. 
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This  procedure  is  only  possible  for  the  nine  contaminants  — 
cadmium,  chromium,  lead,  manganese,  nickel,  nitrogen  oxides, 
nitrous  oxide,,'  sulphur  dioxide  and  PM  —  for  which  monitoring 
data  are  available.-  The  procedure  was  not  used  with  VOC's 
because  the  monitoring  stations  report  total  VOC's  and  we  did  not 
have  a  reasonable  basis  for  speciating  this  total  at  each 
monitoring  station. 

G.4.4   Maps  of  Environmental  Resources  at  Risk 

Having  determined  the  geographic  dispersion  of  a  contaminant  and 
its  concentration  at  each  point,  it  is  necessary  to  relate  that 
information  to  the  resources  at  risk.   These  include  humans, 
crops,  livestock,  forests  and  surface  water.   Maps  of  these 
resources  were  purchased  from  various  sources  and  converted  to 
SPANS  format. 

Population  by  enumeration  area  for  the  1986  census  was  purchased 
from  Statistics  Canada.   The  enumeration  area  population  totals 
were  aggregated  to  the  census  subdivision  level.   The  area  of 
each  census  subdivision  was  calculated  by  SPANS.   This  allowed 
the  population  figure  to  be  converted  to  a  population  density 
(people  per  sq.  km) .   The  population  density  is  assumed  to  be 
uniform  across  the  census  subdivision. 

Agricultural  data  from  the  1986  Census  of  Agriculture  were 
purchased  from  Statistics  Canada.   The  information  is  available 
at  the  consolidated  census  sub-division  level.   Only  the  total 
crop  acreage  by  consolidated  census  sub-division  was  mapped. 

Forestry  data  were  obtained  from  the  National  Forestry  Service  in 
Petawawa.   A  map  of  surface  water  for  Ontario  was  purchased  from 
the  Canada  Lands  Directorate  Service. 
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G.4.5   Analysis  of  Environmental  Resources  at  Risk 

For  most  contaminants,  only  the,  estimated  emissions  due  to 
Regulation  StisVsouirces  are  available.   In  these  cases  where  the 
analysis  can  only  identify  the  environmental  resources  at  risk. 
The  dispersion  pattern  for  the  contaminant  is  superimposed  on  the 
appropriate  environmental  resource(s),  usually  population.   This 
is  done  for  each  scenario.   The  difference  between  the  current 
resource  exposure  and  the  exposure  under  each  scenario  is 
calculated.   This  indicates  the  maximum  change  in  environmental 
resources  at  risk  by  scenario. 

In  each  of  these  cases,  the  estimated  emissions  due  to 
Regulation  308  sources  are  also  compared  to  the  threshold  values 
above  which  observable  effects  occur.   In  virtually  every 
instance  the  current  emissions  are  well  below  the  threshold 
level. 

G.4.6   Analysis  of  Contaminants  for  which  Ambient  Concentrations 
have  been  Estimated 

The  nine  contaminants  —  nickel,  chromium,  lead,  manganese, 
nickel,  nitrogen  oxides,  nitrous  oxide,  sulphur  dioxide  and  PM 
—  for  which  total  ambient  concentrations  are  estimated  are 
analysed  in  a  manner  similar  to  that  outlined  above. 

The  areas  where  the  total  ambient  concentration  exceeds  the 
threshold  value  for  observable  effects  is  plotted  for  the 
existing  situation  and  for  each  scenario  for  implementation  of 
proposed  changes  to  Regulation  308.   The  difference  in  the 
resources  (population,  crops,  forests)  exposed  to  above-threshold 
levels  of  the  contaminant  is  a  measure  of  the  public  benefits  of 
the  proposed  revisions  to  Regulation  308.  * 
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G.4.7   Analysis  of  Cases  with  Exposure-Response  and  Damage 
Functions 

In  the  case  of'  cancer,  and  systemic  health,  visibility  and 
soiling  effects  in  the  case  of  a  few  important  contaminants, 
exposure-response  and  damage  functions  are  available.   These 
functions  calculate  the  benefits  of  reduced  exposure  from  the 
change  in  exposure  and  the  population  exposed. 

For  each  scenario,  the  reduction  in  exposure  from  the  existing 
level  is  calculated  at  each  grid  point.   This  change  is 
multiplied  by  the  population  density  in  the  4  km  square  of  which 
the  grid  point  is  the  centre  and  by  the  area  (16  km^)  to  get  the 
population  and  change  in  exposure.   These  values  are  summed  over 
all  relevant  grid  points  and  the  aggregate  values  are  used  in  the 
exposure-response  and  damage  function  to  calculate  the  estimated 
benefits  of  reduced  emissions. 

Maps  showing  the  dispersion  patterns  by  contaminant  are  provided 
in  Appendix  I. 
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H.l   INTRODUCTION 


An  exposure-response  function  predicts  how  health  and  other 
factors,  such-'as  visibility  and  materials  damage,  change  as  the 
concentration  of  a  contaminant  changes.   Exposure-response 
functions  are  the  most  commonly  available  for  deaths  by  cancer. 
Functions  are  also  available  for  other  health  effects  from 
nitrogen  dioxide,  sulphur  dioxide,  and  suspended  particular 
matter. 

The  exposure-response  functions  are  discussed  below  for  cancer 
and  various  effects  of  nitrogen  dioxide,  sulphur  dioxide,  and 
suspended  particulate  matter.   The  functions  used  are  listed 
under  the  contaminant  in  question  in  Appendix  C.   They  are 
elaborated  here  for  greater  clarity.   For  many  non-cancer 
effects,  upper,  lower,  and  central  values  for  the  function  are 
listed  in  Appendix  C.   We  have  calculated  the  lower  and  central 
value  for  each  function. 


H.2   CANCER  FUNCTIONS 


The  change  in  cancer  deaths  due  to  a  change  in  the  concentration 
of  a  contaminant  is  predicted  using  a  unit  risk  factor.   The 
unit  risk  factors  we  used  are  usually  IRIS  (Integrated  Risk 
Information  System)  values.   Where  IRIS  values  are  not 
available,  we  used  Massachusetts  unit  risk  factors. 

The  change  in  annual  cancer  mortality  from  a  change  in  the 
concentration  of  a  chemical  is  calculated  using  the  following 
formula : 


H  -  2 
ACM  =  URF  X  CONC  X  POP  /  70 

Where:        •; 

ACM   =   Change  in. annual  cancer  mortality  resulting  from  the 

change  in  concentration  of  the  contaminant. 
URF    =   The  IRIS  (or  Massachusetts)  Unit  Risk  Factor. 
CONC  =   Change  in  concentration  of  the  contaminant  measured  on 

an  annual  basis  expressed  in  micrograms  per  cubic 

metre. 
POP   =   The  population  subject  to  the  change  in  contaminant 

concentration. 
70    =   A  factor  converting  lifetime  (estimated  to  be  70 

years)  deaths  to  annual  deaths. 

The  change  in  annual  cancer  mortality  (ACM)  is  calculated  for 
each  contaminant  for  which  a  cancer  unit  risk  factor  is 
available.   The  calculation  is  performed  for  each  of  the 
implementation  scenarios. 

The  calculation  is  performed  as  follows  for  a  given  contaminant 
and  scenario.   Each  of  the  over  15,000  grid  cells  used  in  the 
analysis  has  a  change  in  contaminant  concentration  specific  to 
the  contaminant  and  scenario  (current  concentration  in  micrograms 
per  cubic  meter  as  an  annual  average  less  projected  concentration 
under  the  specified  scenario  in  micrograms  per  cubic  metre  as  an 
annual  average).   The  change  in  contaminant  concentration  (CONC) 
for  a  grid  cell  is  multiplied  by  the  population  (POP)  of  the  grid 
cell.   The  products  of  these  multiplications  are  summed  over  all 
grid  cells.   The  sum  is  then  multiplied  by  the  cancer  unit  risk 
factor  (URF)  and  divided  by  7  0  to  yield  the  change  in  annual 
cancer  mortality. 

Using  arsenic  as  an  example,  the  change  in  concentration 
multiplied  by  the  population  summed  over  all  grid  cells  is  as 
follows  by  scenario: 
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Change  in    Unit  Risk 
Concentration   Factor 
X  Population 


Change  in 
Annual  Cancer 
Mortality* 


Scenarios  "A"  &  "D" 

117, 

,335 

4. 

.3 

X 

10--' 

7.2 

Scenario  "B" 

128, 

,240 

4. 

.3 

X 

10-3 

7.9 

Scenario  "C" 

105, 

,577 

4, 

,3 

X 

10-3 

6.5 

Scenario  "E" 

113, 

,821 

4, 

.3 

X 

10-3 

7.0 

117,335  X  0.0043/70  =  7.2 

The  unit  risk  factor  for  each  carcinogen  for  which  information  is 
available  is  shown  in  Exhibit  H.l. 


H.3   NITROGEN  DIOXIDE  FUNCTIONS 


Exposure-response  functions  for  nitrogen  dioxide  have  been 
developed  for: 

.  eye  irritation; 

.  visibility  effects;   and 

.  materials  damage. 

The  functions  listed  below  are  derived  from  the  material 
presented  in  Appendix  C. 

H.3.1   Eve  Irritation 


The  eye  irritation  damage  function  depends  upon  measurement  of 
average  daily  high  concentration  of  NO2 .   We  only  had  information 
on  annual  average  concentration  and  therefore  were  unable  to 
calculate  changes  eye  irritation  under  the  proposed 
implementation  scenarios  due  to  reduced  emissions  of  NO2 . 
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EXHIBIT  H.l:  '  CANCER  UNIT  RISK  FACTORS 


Contan^inant  Unit  Risk  Factor 

Acrylonitrile  6.8  x  10"^ 

Arsenic  4.3  x  lO'^ 


-3 
-6 


Asbestos  7.6  x  10 

Benzene  8.3  x  10 

1,3  Butadiene  2.8  x  10""^ 


Cadmium  1.8  x  10 


-4 


Carbon  Tetrachloride  1.5  x  10"^ 

Chloroform  2.3  x  10"^ 

Chromium  1.2  x  IC^ 

Epichlorohydrin  1.2  x  10"^ 

Ethylene  Bichloride  2.6  x  10"^ 

Formaldehyde  1.3  x  10"^ 

Methylene  Chloride  4.1  x  10"^ 


Nickel  4.8  X  10 


-4 


Polychlorinated  Biphenols  2.2  x  10"-^ 

Propylene  Oxide  6.67  x  10"'^ 

Styrene  5.7  x  10""^ 

Toluene  Diisocyanate  6.79  x  10"^ 
1,1, 1-Trichloroethane 


(Methyl  Chloroform)  1.6  x  10 


-5 


Trichloroethylene  1.3  x  10"^ 

Vinyl  Chloride  2.6  x  10"^ 


Source:   Appendix  C 


H  -  5 
H.3.2   Visibility  Effects 

We  adjusted  the  function  given  in  Appendix  C  as  follows: 

*'"'■'  •    '. 
Total  ■  "  " 

Value      =     Ex  205  x  In  [l/(l-{0.035  x  (NO2A  "  NO2B)  x 
Per  Year       IOO/NO2A))]  x  POP/HF 

Where: 

E      =   The  mean  1986  rate  of  exchange  (1.3894)  to  convert 

1986  U.S.  to  1986  Canadian  dollars. 
In     =  The  natural  logarithm. 
0.035   =  A  mean  value  for  the  function,  which  varies  according 

to  urban  or  rural  location.   (Taking  a  mean  value 

biases  the  result,  but  the  difference  between  urban  and 

rural  location  is  not  as  great  as  the  upper  and  lower 

bounds  for  this  function.) 
N'-'ZA   ~  "^^^  current  NO2  concentration  for  the  grid  cell. 
NO2B   =  The  NO2  concentration  for  the  same  grid  cell  under  the 

specified  implementation  scenario, 
POP    =   The  population  of  the  grid  cell. 
HF     =  A  factor  (2.8)  to  convert  population  to  the  number  of 

households. 

The  function  given  above  is  the  central  case  estimate.   To  get 
the  lower  bound  estimate,  the  constant  205  is  replaced  by  100. 
In  other  words,  the  lower  bound  is  equal  to  the  central  case 
divided  by  2.05, 

As  was  the  case  for  the  cancer  mortality  equation,  the 
calculation  is  performed  for  each  of  the  over  15,000  grid  cells 
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and  then  summed.   Specifically, 

In  [l/(l-(0.035  X  (NO2A  -  NO2B)  X  IOO/NO2A))]  X  POP/HF 

is  calculated  for  each  grid  cell.   These  values  are  then  summed 
over  all  grid  cells.   And  the  sum  is  multiplied  by  E  x  205  to  get 
the  final  result. 

Note  that  the  output  is  a  dollar  value;   there  is  no  physical 
value  such  as  "number  of  days  of  poor  visibility",  analogous  to 
the  number  of  deaths,  hospital  days,  etc.  in  the  other  functions. 

H.3.3   Materials  Damage 

The  materials  damage  function  for  nitrogen  dioxide  is  expressed 
in  similar  terms.   The  equation  describes  the  annual  materials 
damage  value  per  person.   For  the  central  case,  the  materials 
damage  equation  is: 

Total 

Value     =  E  X  1.51  X  [ (NO2A  -  NO2B)  x  IOO/NO2A]  x  POP 

Per  Year 

The  symbols  have  the  same  meaning  as  in  the  visibility  equation. 
Thus,  materials  damage  is  a  function  of  the  percentage  change  in 
nitrogen  dioxide  concentration  and  population. 

The  function  shown  above  is  the  central  case.   The  lower  bound  is 
calculated  as  0.5  times  the  central  case. 
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H.4   SULPHUR  DIOXIDE  FUNCTIONS 


H.4.1   Mortal-itv 

An  exposure-response  function  is  available  for  deaths 
attributable  to  sulphur  dioxide.   The  function  in  Appendix  C  is 
adjusted  as  follows  to  calculate  absolute  change  in  the  number  of 
annual  deaths. 

Deaths  =  2.672  X  10"^  x  POP  x  (SO2A-SO2B) 

Where : 

^"^ZA   ~   '^^^  ^^2  concentration  for  a  grid  cell  under  existing 

regulations. 
SO2B   =  The  SO2  concentration  for  a  grid  cell  under  a  specific 

implementation  scenario. 
POP    =  The  population  of  a  grid  cell. 

The  function  presented  above  is  the  central  case  estimate  for 
reduction  in  annual  mortality  due  to  reduced  concentrations  of 
sulphur  dioxide.   The  low  case  is  calculated  by  applying  the 
lower  bound  estimate  of  the  value  of  a  statistical  life  to  the 
calculated  reduction  in  mortality. 

H.4. 2   Hospital  Days  for  Respiratory  Conditions 

The  function  of  hospital  days  for  respiratory  conditions  (HDRC) 
is: 

HDRC   =    4.0  X  10"'*  x  (SO2A  -  SO2B)  ^    POP 

The  variables  S02a»  SO2B  and  POP  have  the  same  meaning  as  in  the 
mortality  equation. 
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The  lower  bound  estimate  is  obtained  by  replacing  4.0  x  lO"'*  in 
the  central  case  equation  presented  above  with  1.5  x  10""*.   In 
other  words 'the  lower  bound  is  0.375  times  the  central  case. 

H.4.3   Hospital  Admissions  for  Respiratory  Disease 

The  function  for  hospital  admissions  for  respiratory  disease 
(HARD)  is: 

HARD    =   4.065  X  10"^  x  (SO2A  "  SO2B)  X  POP 

The  variables  SO2A/  SO2B  arid  POP  have  the  same  meaning  as  in  the 
mortality  equation. 

The  lower  bound  estimate  is  calculated  by  replacing  4.065  x  10"^ 
in  the  central  case  equation  presented  above  with  0.63  x  10"^. 
Alternately,  the  lower  bound  is  0.155  times  the  central  case. 


H.5   PARTICULATE  MATTER  FUNCTIONS 


Six  exposure  response  functions  exist  for  PM.   They  relate  to 

.  death; 

.  emergency  room  visits; 

.  restricted  activity  days; 

.  chronic  obstructive  pulmonary  disease; 

.  visibility;   and 

.  materials  damage. 
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The  functions  below  are  derived  from  those  in  Appendix  C.   In 
each  case: 

PMj^    =  The'original  PM  concentration  for  a  grid  cell. 
PMg    =  The  PM  concentration  for  a  grid  cell  under  the 

scenario  considered. 
E      =  The  mean  1986  rate  of  exchange  (1.3894)  to  convert 

1986  U.S.  to  1986  Canadian  dollars. 
HF     =  A  factor  (2.8)  to  convert  the  population  of  a  grid 

cell  to  the  number  of  households. 
POP    =  The  population  of  a  grid  cell. 

H.5. 1   Death 

Deaths  =   2.21  X  10"^  x  POP  x  (PM^^-PMg) 

This  equation  gives  the  change  in  annual  deaths  due  to  a  change 
in  particulate  matter  concentration.   The  lower  bound  and  central 
value  are  calculated  by  applying  the  corresponding  values  of  a 
statistical  life  to  the  reduction  in  annual  mortality  estimated 
using  the  above  equation. 

H.5.2   Emergency  Room  Visits 

ERV  =  1.3  X  10""^  X  POP  X  (PM^-PMb) 

The  change  in  the  number  of  emergency  room  visits  per  year,  ERV, 

is  related  to  the  change  in  particulate  matter  concentration. 

The  lower  bound  estimate  of  emergency  room  visits  is  calculated 

by  replacing  the  central  case  coefficient  of  1.3  x  lO"**  by 
1.2  X  10~5. 
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H.5.3   Restricted  Activity  Days 

RAD  =   2.736  X  10"^  X  POP  X  (PMj^-PMg)  -  Emergency  Room 

Visits 

The  change  in  the  number  of  restricted  activity  days,  RAD,  is 
reduced  by  the  change  in  the  number  of  emergency  room  visits 
(ERV)  to  avoid  double  counting.   The  lower  bound  estimate  of 
restricted  activity  days  is  calculated  by  replacing  the  central 
case  coefficient  of  2.736  x  10"^  by  1.938  x  lO"^.   in  other 
words,  the  lower  bound  is  equal  to  0.70833  times  the  central 
case. 

H.5.4   Chronic  Obstructive  Pulmonary  Disease 

The  function  in  Appendix  C  deals  with  changes  in  risk  of 
pulmonary  disease  in  adults  and  children.   It  does  not  predict 
actual  change  in,  or  values  relating  to,  pulmonary  disease.   It 
was  therefore  impractical  to  use  it  in  the  analysis. 

H.5.5   Visibility  Effects 

Total 

Value  =Ex206xln    (PM^/PMb)    x    POP/HF 

Per  Year 

The  lower  bound  estimate  is  calculated  as  0.5  times  the  central 
case. 

As  with  visibility  effects  from  nitrogen  dioxide,  this  function 
returns  an  annual  dollar  value  for  improved  visibility. 

H.5.6   Materials  Damage 

Total 

Value      =  E  X  7.4  X  (PM^^-PHg)  x  POP/HF 

Per  Year 
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This  function. also  returns  an  annual  dollar  value  of  reduced 
materials  damage. 

To  get  the  lower  bound  estimate,  the  central  case  coefficient  of 
7.4  is  replaced  by  0.52.   In  other  words,  the  lower  bound  is 
0.07027  times  the  central  case. 
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EXHIBIT  1.2:   LOCATION  OF  THE  LINE  ALONG  WHICH  THE  CONCENTRATION 
PROFILE  IS  CALCULATED 
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PROFILE  OF  ARSENIC  CONCENTRATION  DUE  TO  REGULATION 
3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN  ONTARIO 
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EXHIBIT  1.6:    PROFILE  OF  ASBESTOS  CONCENTRATION  DUE  TO  REGULATION 
3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN  ONTARIO 
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EXHIBIT  1.8:    PROFILE  OF  BENZENE  CONCENTRATION  DUE  TO  REGULATION 
3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN  ONTARIO 
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EXHIBIT  I. 10:   PROFILE  OF  1,3  BUTADIENE  CONCENTRATION  DUE  TO 

REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 
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EXHIBIT  1.12:   PROFILE  OF  2-BUTANONE  CONCENTRATION  DUE  TO 

REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 
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EXHIBIT  1.14 


PROFILE  OF  BUTYL  ACRYLATE  CONCENTRATION  DUE  TO 
REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 
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EXHIBIT  1.16:   PROFILE  OF  CADMIUM  CONCENTRATION  DUE  TO  REGULATION 
3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN  ONTARIO 
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EXHIBIT  1.18: 


PROFILE  OF  CARBON  TETRACHLORIDE  CONCENTRATION  DUE 
TO  REGULATION  308  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 
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EXHIBIT  1.20:   PROFILE  OF  CDD  CONCENTRATION  DUE  TO  REGULATION  308 
SOURCE  EMISSIONS  ACROSS  SOUTHERN  ONTARIO 
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EXHIBIT  1.22 


PROFILE  OF  CDF  CONCENTRATION  DUE  TO  REGULATION  3  08 
SOURCE  EMISSIONS  ACROSS  SOUTHERN  ONTARIO 
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EXHIBIT  1.24:   PROFILE  OF  CHLORINE  CONCE;:  -.ATI 
308  SOURCE  EMISSIONS  ACROr '  301 
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EXHIBIT  1.26; 


PROFILE  OF  CHLORINE  DIOXIDE  CONCENTRATION  DUE  TO 
REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 
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PROFILE  OF  CHLOROFORM  CONCENTRATION  DUE  TO 
REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 
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EXHIBIT  1.30:   PROFILE  OF  CHROMIUM  CONCENTRATION  DUE  TO  REGULATION 
3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN  ONTARIO 
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EXHIBIT  1.32: 


PROFILE  OF  EPICHLOROHYDRIN  CONCENTRATION  DUE  TO 
REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 
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EXHIBIT  1.34:   PROFILE  OF  ETHYL  BENZENE  CONCENTRATION  DUE  TO 

REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
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EXHIBIT  1.36:   PROFILE  OF  ETHYLENE  CONCENTRATION  DUE  TO  REGULATION 
3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN  ONTARIO 
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Scenario  C 
Scenario  E 


1-37 


c 


cr 

in 

C 

o 
in 

i. 


in  IS  s  s  i/i 
m  ®  s  ■« 

I      I      I     n     01 


-  VT 
X  z 
i;  c 


1-38 


EXHIBIT  I.: 


ETHYLENE  DIBROMIDE  CONCENTRATION  DUE  TO 
3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 


Ethylene  Dibromide 


0.C4 


0.03 


o.c: 


0.01 


0- 


. 

u  ._ 

-80        -79        -78        -77        -76        -75        -74 
Longitude 

West  -  East 


Legend 


Exis: 

Sce^ 
Scene 
Scene- 
Scena- 
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EXHIBIT  1.40:   PROFILE  OF  ETHYLENE  BICHLORIDE  CONCENTRATION  DUE  TO 
REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Ethylene  Dichloride 


Concentration 
(ug/m3 


0.005 


0  1 .-r^— ^-vt^--^i — i-^=^ ^1 i-T 1 r-^5w — ,-a^ 

-84       -83        -82        -81        -80        -79        -78        -77        -76        -75        -74 

Longitude 
West  -  East 
Health  Threshold  •  80  ug/m3 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.42 


1-42 


PROFILE  OF  BUTYL  CELLOSOLVE  CONCENTRATION  DUE  TO 
REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Butyl  Cellosolve 


0.14 


Concentration 
(ug/m3 


0.02 


-84       -83       -82 


Health  Threshold  •  25  ug/m3 


-80        -79        -78 

Longitude 

West  -  East 


77 


-76 


-75        -74 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.44:   PROFILE  OF  CELLOSOLVE  CONCENTRATION  DUE  TO 

REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


30 


25 


20 


15 


10 


Concentration 
(ug/m3)  (1E-5) 


l\ 


Cellosolve 


-84        -83        -82        -81 


Health  Threshold  •  55.4  ug/m3 


-80        -79        -78        -77        -76 
Longitude 
West  -  East 


-75        -74 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.46:   PROFILE  OF  CELLOSOLVE  ACETATE  CONCENTRATION  DUE  TO 
REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Cellosolve  Acetate 


Concentration 
(ug/m3)  (1E-5) 


-84        -83        -82         -81        -80        -79        -78        -77        -76        -75        -74 

Longitude 
West  -  East 
Health  Threshold  •  IS.O  ug/ni3 


Legend 


Existing 

Scenarios  a  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.48:   PROFILE  OF  ETHYLENE  OXIDE  CONCENTRATION  DUE  TO 
REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Ethylene  Oxide 


30 


25 


20 


15 


lOr 


Concentration 
(ug/m3)  (1E-4) 


1 

\ 

.     L  ___ 

cT:t. 

.  «o.,  < 

J  A 

-84        -83        -82         -81        -80        -79        -78        -77        -76        -75        -74 

Longitude 
West  -  East 
Health  Threshold  •  1.0  ug/m3 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 


1-49 


"^ 

© 

\ 

1-50 


EXHIBIT  1.50:   PROFILE  OF  FLUORIDES  CONCENTRATION  DUE  TO 

REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Fluorides 


Concentration 
(ug/m3) 


0.06 


0.05 


0.04 


0.03  - 


0.02  - 


0.01 


0  I ■■..■■-■    I  ■  .pi    l-'l  ,  ,  ,ri'- ^fi^Tt   I  ipA.. 

-84        -83        -82        -81        -80        -79        -78        -77        -76        -75        -74 

Longitude 
West  -  East 
Health  Threshold  •  0.344  ug/m3 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.52:   PROFILE  OF  FORMALDEHYDE  CONCENTRATION  DUE  TO 

REGULATION  308  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Formaldehyde 


0.12 


0.1 


0.08 


0.06 


0.04 


0.02 


Concentration 
(ug/m3) 


- — ^ttA ■ 


-84        -83        -82        -81        -80        -79        -78        -77        -76        -75        -74 

Longitude 
West  -  East 
Health  Threshold  •  5.42  ug/m3 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.54:   PROFILE  OF  HYDROGEN  CHLORIDE  CONCENTRATION  DUE  TO 
REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Hydrogen  Chloride 


Concentration 
(ug/m3) 


0.4 


0.3 


0.2 


0.1  - 


-84        -83        -82        -81        -80        -79        -78        -77        -76        -75        -74 

Longitude 
West  -  East 

Health  Threshold  •  8.0  ug/m3 


Legend 


>\ ^^ .   

Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.56:   PROFILE  OF  LEAD  CONCENTRATION  DUE  TO  REGULATION  308 
SOURCE  EMISSIONS  ACROSS  SOUTHERN  ONTARIO 


Lead 


Concentration 
(ug/m3) 


-84       -83       -82        -81       -80       -79       -78       -77       -76       -75       -74 

Longitude 
West  -  East 
Health  Threshold  •  1.0  ug/m3 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.58:   PROFILE  OF  MANGANESE  COMPOUNDS  CONCENTRATION  DUE  TO 
REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Manganese  Compounds 


0.12 


Concentration 
(ug/m3) 


Health  Threshold  •  2.0  ug/m3 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.60:   PROFILE  OF  MERCURY  CONCENTRATION  DUE  TO  REGULATION 
308  SOURCE  EMISSIONS  ACROSS  SOUTHERN  ONTARIO 


Mercury 


0.005 


0.004 


0.003 


0.002 


0.001 


Concentration 
(ug/m3) 


i 

_           .          y\    y^^J^ 

\ ^— ^ 1., 

-84       -83        -82        -81        -80        -79        -78        -77        -76        -75        -74 

Longitude 
West  -  East 
Health  Threshold  •  C.4  uB/m3 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.62:   PROFILE  OF  METHYL  ISOBUTYL  KETONE  CONCENTRATION  DUE 
TO  REGULATION  308  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Methyl  Isobutyl  Ketone 


Concentration 
(ug/m3) 


-83        -82 


-81 


-80        -79        -78 

Longitude 

West  -  East 


-77 


-76        -75 


-74 


Health  Threshold  ■  240  ug/m3 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.64:   PROFILE  OF  METHYLENE  CHLORIDE  CONCENTRATION  DUE  TO 
REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Methylene  Chloride 


50 


40 


30 


20 


Concentration 
(ug/m3)  (1E-5) 


10 


r7^ ^^n^::^, 


-84        -83        -82         -81        -80        -79        -78        -77        -76        -75        -74 

Longitude 
West  -  East 
Health  Threshold  •  583  ug/m3 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  3 
Scenario  C 
Scenario  E 
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EXHIBIT  1.66:   PROFILE  OF  NAPHTHALENE  CONCENTRATION  DUE  TO 

REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Concentration 
{ug/m3)  (1E-5) 


Naphthalene 


-84        -83        -82        -81 


Health  Threshold  •  2.5  ug/m3 


-80        -79        -78 

Longitude 

West  -  East 


-77        -76        -75 


-74 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.68:   PROFILE  OF  NICKEL  CONCENTRATION  DUE  TO  REGULATION 
3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN  ONTARIO 


Nickel 


0.025 


0.02  — 


0.015  — 


0.01  - 


0.005 


Concentration 
(ug/m3) 


-84       -83       -82        -81        -80       -79       -78       -77       -76       -75       -74 


Longitude 
West  -  East 


Health  Threshold  •  0.4  ug/m3 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 


1-69 


1-70 


EXHIBIT  1.70; 


PROFILE  OF  NITROGEN  OXIDES  CONCENTRATION  DUE  TO 
REGULATION  308  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Nitrogen  Oxides 


2.5 


Concentration 
(ug/m3 


Health  Threshold  •  40  ug/rn3 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.72:   PROFILE  OF  PERCHLOROETHYLENE  CONCENTRATION  DUE  TO 
REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Perchloroethylene 


20 


15 


10 


Concentration 
(ug/m3)  (1E-4) 


1       .^_  ^ ^ - 

.:ii.„_  _ 

f\ 

\     A 

-84        -83        -82         -81        -80        -79        -78        -77        -76        -75        -74 

Longitude 
West  -  East 

Health  Threshold  •  800  ug/m3 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.74:   PROFILE  OF  PHENOL  CONCENTRATION  DUE  TO  REGULATION 
3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN  ONTARIO 


Phenol 


0.4 


Concentration 
(ug/m3) 


Health  Threshold  •  20  ug/m3 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.76:   PROFILE  OF  POLY CHLORINATED  BIPHENYLS  CONCENTRATION 
DUE  TO  REGULATION  3  08  SOURCE  EMISSIONS  ACROSS 
SOUTHERN  ONTARIO 


Polychlorinated  Biphenyls 


30 


25 


20 


15 


10 


Concentration 
(ug/m3)  {1E-5) 


-84        -83        -82        -81 


-80        -79        -78 

Longitude 

West  -  East 


-77        -76 


-75 


-74 


Health  Threshold  •  0.03  ug/m3 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1-78:  PROFILE  OF  PROPIONALDEHYDE  CONCENTRATION  DUE  TO 
REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Concentration 
(ug/nn3)  (lE-5) 


Propionaldehyde 


-84        -83        -82        -81 


Health  Threshold  •  0.5  ug/m3 


-80        -79        -78        -77 
Longitude 
West  -  East 


-76         -75         -74 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.80:   PROFILE  OF  PROPIONIC  ACID  CONCENTRATION  DUE  TO 
REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Concentration 
(ug/m3)  (1E-5 


Propionic  Acid 


-84        -83        -82 


-81 


-80        -79        -78        -77 
Longitude 
West  -  East 


74 


Health  Threshold  •  6.67  ug/m3 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.82:  PROFILE  OF  PROPYLENE  OXIDE  CONCENTRATION  DUE  TO 
REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Propylene  Oxide 


Concentration 
(ug/m3)  (1E-4) 


-84        -83        -82 


-81        -80        -79        -78 
Longitude 
West  -  East 


-77 


-76 


Health  Threshold  •  900  ug/m3 


Legend 


Existing 

Scenarios  A  a  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.84:   PROFILE  OF  SELENIUM  CONCENTRATION  DUE  TO  REGULATION 
3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN  ONTARIO 


Selenium 


Concentration 
(ug/m3)  {1E-4) 


Health  Threshold  ■  2.0  ug/m3 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.86:   PROFILE  OF  SILICA  CONCENTRATION  DUE  TO  REGULATION 
3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN  ONTARIO 


Silica 


1.4 


Concentration 
(ug/m3) 


-84        -83        -82         -81         -80        -79        -78        -77        -76        -75 

Longitude 
West  -  East 
Health  Threshold  ■  1.0  ug/m3 


-74 


Legend 


Existing 

Scenarios  a  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.88: 


PROFILE  OF  STYRENE  CONCENTRATION  DUE  TO  REGULATION 
3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN  ONTARIO 


Styrene 


Concentration 
(ug/m3) 


\\ 

V 

r^         .                  /         .?::■.- 

0.8 


0.6 


0.4 


0.2  - 


-84        -83        -82         -81         -80        -79        -78        -77        -76        -75        -74 

Longitude 
West  -  East 
Health  Threshold  •  80  ug/m3 


Legend 


Existing 

Scenarios  A  a  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.90:  PROFILE  OF  SULPHUR  DIOXIDE  CONCENTRATION  DUE  TO 
REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Sulphur  Dioxide 


Concentration 
(ug/m3) 


-84        -83        -82   ■     -81        -80        -79        -78        -77        -76        -75        -74 


Longitude 
West  -  East 


Health  Threshold  -  55  ug/m3 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.92:   PROFILE  OF  TOTAL  PARTICULATE  CONCENTRATION  DUE  TO 
REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 

ONTARIO 


Total  Particulate 


Concentration 
(ug/m3) 


-80        -79        -75 

Longitude 

West  -  East 


74 


i:.:\b  Threclcld  ■  24  ug. 


Legend 


rxisting 
Scenarios  - 
Scenario  5 
Scenario  C 
Scenario  E 
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EXHIBIT  1.94:   PROFILE  OF  TIN  CONCENTRATION  DUE  TO  REGULATION  308 
SOURCE  EMISSIONS  ACROSS  SOUTHERN  ONTARIO 


Tin 


Concentration 
(ug/m3) 


0.004 


0.003 


0.002 


0.001 


-84       -83       -82        -81        -80       -79       -78       -77       -76       -75       -74 

Longitude 
West  -  East 

Health  Threshold  ■  2.0  ug/m3 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.96:   PROFILE  OF  TOLUENE  CONCENTRATION  DUE  TO  REGULATION 
3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN  ONTARIO 


Toluene 


Concentration 
(ug/m3) 


■84        -83        -82         -81 


Health  Threshold  •  400  ug/m3 


-80        -79        -78 

Longitude 
West  -  East 


-77         -76         -75 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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PROFILE  OF  TOLUENE  DIISOCYANATE  CONCENTRATION  DUE 
TO  REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Toluene  Diisocyanate 


Concentration 
(ug/m3)  (1E-5 


-84 


-83        -82 


n         -80        -79        -78         -77 
Longitude 
West  -  East 


-74 


Health  Threshold  •  C."  _• 


Legend 


Existing 

Scena''cs  A  3  D 
Scenario  B 
Scenario  C 
Scenario  £ 
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EXHIBIT  I. 100:  PROFILE  OF  1, 1, 1-TRICHLOROETHANE  CONCENTRATION  DUE 
TO  REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


1,1,1-Trichloroethane 


Concentration 
(ug/m3) 


0.004 


0.003  - 


0.002  - 


0.001  - 


-84       -83        -82        -81        -80        -79        -78        -77        -76        -75        -74 

Longitude 
West  -  East 

Health  Threshold  •  23000  ug/m3 


Legend 
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Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.104:  PROFILE  OF  TRICHLOROETHYLENE  CONCENTRATION  DUE  TO 
REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Trichloroethylene 


35 
30 
25 
20 
15 
10 
5 


Concentration 

(ug/m3)  (1E-4) 


^  _ J_.._ _ ,  _ . 

1 

^ 



.  ,--. ._ : 

-  ^.M-^-Uy-)r^  .     ,     ,  ^ 

-84        -83        -82        -81 


Health  Threshold  ■  5600  ug/m3 


-80        -79        -78 

Longitude 
West  -  East 


-77 


-76        -75        -74 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.106:  PROFILE  OF  VINYL  CHLORIDE  CONCENTRATION  DUE  TO 
REGULATION  3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN 
ONTARIO 


Vinyl  Chloride 


0.035 


Concentration 
(ug/m3) 


-80       -79       -7{ 

Longitude 

West  -  East 


75 


-74 


Health  Threshold  •  56  ug/m3 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 
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EXHIBIT  1.108:  PROFILE  OF  XYLENES  CONCENTRATION  DUE  TO  REGULATION 
3  08  SOURCE  EMISSIONS  ACROSS  SOUTHERN  ONTARIO 


Xylenes 


Concentration 
(ug/m3)  (1E-5) 


-84        -83        -82        -81        -80        -79        -78        -77        -76        -75 

Longitude 
West  -  East 

Health  Threshold  ■  460  ug/m3 


-74 


Legend 


Existing 

Scenarios  A  &  D 
Scenario  B 
Scenario  C 
Scenario  E 


